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ABSTRACT 

The Wilkinson Microwave Anisotropy Probe ( WMAP) has mapped the entire 
sky in five frequency bands between 23 and 94 GHz with polarization sensitive 
radiometers. We present three-year full-sky maps of the polarization and analyze 
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them for foreground emission and cosmological implications. These observations 
open up a new window for understanding how the universe began and help set a 
foundation for future observations. 

WMAP observes significant levels of polarized foreground emission due to 
both Galactic synchrotron radiation and thermal dust emission. Synchrotron 
radiation is the dominant signal at £ < 50 and v < 40 GHz, while thermal dust 
emission is evident at 94 GHz. The least contaminated channel is at 61 GHz. We 
present a model of polarized foreground emission that captures the large angular 
scale characteristics of the microwave sky. 

After applying a Galactic mask that cuts 25.7% of the sky, we show that the 
high Galactic latitude rms polarized foreground emission, averaged over £ = 4 — 6, 
ranges from as 5 fiK at 22 GHz to < 0.6 fiK at 61 GHz. By comparison, the 
levels of intrinsic CMB polarization for a ACDM model with an optical depth of 
t = 0.09 and assumed tensor to scalar ratio r = 0.3 are ~ 0.3 fiK for E-mode 
polarization and « 0.1 fiK for B-mode polarization. To analyze the maps for 
CMB polarization at £ < 16, we subtract a model of the foreground emission 
that is based primarily on a scaling WMAP's 23 GHz map. 

In the foreground corrected maps, we detect £(£ + l)C^ 2 _ 6> /27r = 0.086 ± 
0.029 (fiK) 2 . This is interpreted as the result of rescattering of the CMB by 
free electrons released during reionization at z r = ll.O^J for a model with 
instantaneous reionization. By computing the likelihood of just the EE data as 
a function of r we find r = 0.10 ± 0.03. When the same EE data are used in the 
full six parameter fit to all WMAP data (TT, TE, EE), we find r = 0.09 ± 0.03. 
Marginalization over the foreground subtraction affects this value by 5t < 0.01. 

We see no evidence for B-modes, limiting them to £(£ + l)Cj?J 3 <2 _ 6> /2n = 
—0.04 ± 0.03 (/iK) 2 . We perform a template fit to the E-mode and B-mode data 
with an approximate model for the tensor scalar ratio. We find that the limit 
from the polarization signals alone is r < 2.2 (95% CL) where r is evaluated 
at k — 0.002 Mpc -1 . This corresponds to a limit on the cosmic density of 
gravitational waves of Qcwh 2 < 5 x 10~ 12 . From the full WMAP analysis, we 
find r < 0.55 (95% CL) corresponding to a limit of VL G wh 2 < 1 x 10~ 12 (95% 
CL). The limit on r is approaching the upper bound of predictions for some of 
the simplest models of inflation, r ~ 0.3. 

Subject headings: cosmic microwave background, polarization, cosmology: obser- 
vations 
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Introduction 



The temperature anisotropy in the cosmic microwave background is well established as 
a powerful constraint on theories of the early universe. A related observable, the polarization 
anisotropy of the CMB, gives us a new window into the physical conditions of that era. At 
large angular scales the polarization has the potential to be a direct probe of the universe 
at an age of 10~ 35 s as well as to inform us about the ionization history of the universe. 
This paper reports on the direct detection of CMB polarization at large angular scales and 
helps set a foundation f or future obse r vation s. It is one of four related papers on the three- 
year WMAP analysi s: iJarosik et al.l (120061 ) report on systematic errors an d mapmaking, 



Hinshaw et al.l (120061 ) on the temperature anisotropy and basic results, and ISpergel et al. 
on the parameter estimation and cosmological significance. 



T he po larization of the CMB was predicted soon after the discovery of the CMB 



(IReesI I19681 ). Since then, considerable advances ha ve been made on both theoretical and 



observational fronts. The theoretical development (Basko Sz Polnarev 198ol; Kaiser 



Bond & Efstathioul ll984J: Polnarev 



Harari fc Zaldarriaga 



199 



3; 



1985 



Frewin et al 



Bond fc Efstathiou 



1994 



198^ 



Coulson et al. 



199- 



Crittenden et al 



Crittenden et al 



19831 : 
19931: 



1995; 



Ng fc Nell995l:IZaldarriaga fc Hararll995l : lKosowskv!ll996l : ISeliakill997l : IZaldarriaga fc Seliak 



1997 



Kamionkowski et al.l 119971 ) has evolved to wh ere there are precise predictions and a 



common language to describe the polarization signal. iHu &: White! (119971 ) give a pedagogical 
overview. 



The first limits on the polarization were placed by 


Penzias & Wilson ( 196. 


5), fc 


)llowed 


bv Caderni et al. (1978): Nanos (1979 


); Lul 


oin &; Smoot 


( 


1979, 


1981): Lubin et 


al.| 


(1983): 


Wollack et al. (1993^: Netterfield et al. 


(1997 


): Sironi et al. 


(1997) 


: Torbet et al. ( 


1999 


); Keating et al. 


(2001) and Hedman et all (20021). In 2002. the DASI team announced a detection of CMB 



polarization at sub- degree angular scales based on 9 m onths of data from a 13 element 30 GHz 
interferometer (IKovac et al.l 120021 : iLeitch et al.l 120021 ) . The signal level was consistent with 
that expected from m easurements of the temperature spectrum. The DASI results were con- 
firme d and extended ( Leitch et al.l 120051 ) al most contemporaneo usly with the release of the 
CBI dReadhead et al.ll2004h and CAPMAP faarkats et al.ll2005h resul ts. More recently, th e 
Boomerang team has released its measurement of CMB polarization (IMontroy et al.l 120051 ). 
All of these measurements were made at small angular scales { I > 100). Of the exper iments 
that measure the polarization, the DASI, CBI, and Boomerang (jPiacentini et al.ll2005l ) teams 
also report detections of the temperature-polarization cross correlation. 

The CMB polarization probes the evolution of the decoupling and reionization epochs. 
The polarization signal is generated by Thompson scattering of a local quadrupolar radiation 
pattern by free electrons. The scattering of the same quadrupolar pattern in a direction 
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perpendicular to the line of sight to the observer has the effect of isotropizing the quadrupolar 
radiation field. The net polarization results from a competition between these two effects. We 
estimate the magnitude of the signal following Basko and Polnarev (1980). By integrating the 
Boltzmann equation for the photon distribution they show that the ratio of the polarization 
anisotropy (E rms ) to the temperature (T rms ) signal in a flat cosmology is given by 

Trms / °°[6e-(*') +e-^')]Vl^ i dz'' 1 ' 

where r(z) = cot f n e (z')dz'(dt/dz r ) is the optical depth. Here, ot is the Thompson cross 
section, c is the speed of light, and n e is the free electron density. The difference in brackets 
in the numerator sets the range in z over which polarization is generated. For example, if 
the decoupling epoch entailed an instantaneous transition from an extremely high optical 
depth (r >> 1) to transparency (r = 0), there would be no polarization signal. 

To estimate the polarization fraction we comput e the optical depth using ordinar y 



atomic physics and the thermal history of the universe (lPeebleslll968l ; IZeldovich et al.lll969l ). 
The result is shown in Figure 1. From inserting r(z) in EquationCD, we find that the expected 
level of polarization anisotropy is ~ 5% (in E rms /T rms ) of the anisotropy. 

The polarization producing quadrupol e is generated by different mechanisms at differ- 



ent epochs. Near decoupling at Zd = 1088 (IPage et al.ll2003bl ; ISpergel et al.l 120031 ). velocity 



gradients in the flow of the primordial plasma give rise to the quadrupole. More specifically, 
in the rest frame of an electron in such a flow, the radiation background has a quadrupolar 
pattern proportional to the velocity gradient, Vv, and the mean free path between scatter- 
ings, A. Just before decoupling, z > Zd, the photons are tightly coupled to the electrons 
and A is small. Thus, the polarization is small. As decoupling proceeds A increases and the 
quadrupole magnitude increases. The process is cut off at lower reds hift because the opti 



cal de pth drops so rapidly. In the context of inflationary cosmology, lHarari &: Zaldarriaga 



(119931 ) show that in Fourier space the polarization signal is cx kvA where k is the wavevector 



and A « A is the width of the last scattering surface. 

After decoupling there are no free electrons to scatter the CMB until the first generation 
of stars ignite and reionize the universe at z r . The free electrons then scatter the intrinsic 
CMB quadrupole, C2(z r ), and produce a polarized signal oc C2(z r ) 1 ' 2 r(z r ). As this process 
occurs well after decoupling, the effects of the scattering are manifest at comparatively lower 
values of I. We expect the maximum value of the signal to be at £ max ~ Tt/Q H (z r ) where 
9h{z t ) is the current angular size of the horizon at reionization. For 6 < z < 30 a simple fit 
gives 9h(z) = 4.8/2; - 7 , so that for z r = 12, £ max « 4. Thus, the signature of reionization in 
polarization is cleanly separable from the signature of decoupling. In the first data release the 
WMAP team published a measurement of the temperature-polarization (TE) cross spectrum 
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for 2 < £ < 4 50 (IBennett et al.l l2003bl ; iKogut et al.l 120031 ) with distinctive anti-peak and 
peak structure (IPage et al.ll2003bl ). The £ > 16 part of the spectrum was consistent with the 
prediction from the temperature power spectrum, while the £ < 16 part showed an excess 
that was interpreted as reionization at 11 < z r < 30 (95% CL). 

This paper builds on and extends these results. Not only are there three times as much 
data, but the analysis has improved significantly: 1) The polarization mapmaking pipeline 
now self-consistently includes almost all known effects and correlations due to instrumenta l 
systematics, gain and offset drifts, unequal weighting, and masking (j Jarosik et al.l 120061 ) . 
For example, the noise matrix is no longer taken to be diagonal in pixel space, leading 
to new estimates of the uncertainties. 2) The polarization power spectrum estimate now 
consistently includes th e temperature, E an d B modes (defined below), and the coupling 
between them (see also iHinshaw et al.l 120061 ). 3) The polarized foreground emission is now 
modeled and subtracted in pixel space ( §4.31) . Potential residual contamination is examined 
£ by £ as a function of frequency. In addition to enabling the production of full sky maps 
of the polarization and their power spectra, the combination of these three improvements 
has led to a new measure of the £ < 16 TE and EE spectra, and therefore a new evaluation 
of the optical depth based primarily on EE. The rest of the paper is organized as follows: 
we discuss the measurement in £J2] and consider systematic errors and maps in §3j In §H we 
discuss foreground emission. We then consider, in £J5] and $61 the polarization power spectra 
and their cosmological implications. We conclude in £0 



The Measurement 



WMAP measures the difference in intensity between t wo beams separated by ~ 140° in 



five frequency bands centere d on 23, 33, 41, 61, and 94 GHz (IBennett et al.ll2003bl : IPage et al. 



2003bl : I Jarosik et al.l l2003ai). These are called K, Ka, Q, V, and W bands respectively. 
Corrugated feeds (IBarnes et al.l 120021 ) couple radiation from back-to-back telescopes to the 
differential radiometers. Each feed supports two orthogonal polarizations aligned so that 
the unit vectors along the direction of maximum electric field for an A-side feed f ollow 



(x s ,y s ,z s ) ~ (±1, — sin 20°, — cos20°)/\/2 in spacecraft coordinates (IPage et al.ll2003bl ). For 
a B-side feed, the directions are (x s , y S} z s ) ~ (±1, sin 20°, — cos 20°)/ v / 2. The z s axis points 
toward the Sun along the spacecraft spi n axis; the y x — z x pl ane bisects the telescopes and 
is perpendicular to the radiator panels (IBennett et al.ll2003bl . Figure 2) (jPage et al.ll2003bl . 
Figure 1). The angle between the spacecraft spin axis and the optical axes is ~ 70°. Thus 
the two polarization axes on one side are oriented roughly ±45° with respect to the spin 
axis. 



- 6 - 



The polarization maps are derived from the difference of two differential measurements 



( Jarosik et al. 


2006; 


Kogut et al.l 


2003; 



Hinshaw et al.l l2003bl ) . One half of one differencing 



assembly (DA) (I Jarosik et al.ll2003al ) measures the difference between two similarly oriented 
polarizations, AT 1; from one feed on the A side and one feed on the B side (e.g., W41: 
polarization 1 of the 4th W-band DA corresponding to x s = +1 in both expressions above). 
The other half of the DA measures the difference between the other polarizations in the same 
pair of feeds, AT 2 (e.g., W42: polarization 2 of the 4th W-band DA corresponding to x s = — 1 
in both the expression above). The polarization signal is proportional to AT X — AT 2 . In other 
words, WMAP measures a double difference in polarized intensity, not the int ensity of the 



difference of electric fields as with interferomete rs and correlation receivers (e.g.. lLeitch et al. 



20021 ; Keating et al.ll200ll ; iHedman et al.ll2002[ ). 



With these conventions , the total intensity and polarization signals as measured at the 
output of the detectors are (IKogut et al.ll2003l . Eq. 3&4): 



AT/ = -(ATi + AT 2 ) 



I in a) - I {jib) 



AT P ee -(AT\-AT 2 ) 

= Q(h A ) cos 27a + U (h A ) sin 27,4 
-Q{h B ) cos 2^ B - U (h B ) sin 27^. 



(2) 

(3) 
(4) 



where tia and n B are the unit vectors for the A and B sides; /, Q, and U are the Stokes 
parameter^, and 7 is the angle between the polarization direction of t he electric field and 



the Galactic meridian (IKogut et al.l 20031) . In the mapmaking algorithm (IWright et al.lll996 
Hinshaw et al.l l2003bl ; I Jarosik et al.l 120061 ) . /, Q, and U maps of the sky are produced from 
the time-ordered differential measurements, AT/ and ATp. From these, we form maps of 
polarization intensity, P = \/ Q 2 + U 2 , and direction, 7 = | tan' 1 (U /Q). This convention 
has 7 positive for No rth through West and follows the convention in IZaldarriaga Seljak 
Jl997h and HEALPix iGorski et alJl998h . however, it differs from the standard astronomical 
position angle (PA) which has 7pa = | tan _1 (— U/Q) with •ypA positive for North through 
East. The choice of convention does not affect the plots. 



For linear polarization in a given pixel, the Q and U quantities are rel ated to the x 
and y components of the electric field, E x ,E y , through the coherency matrix (IBorn fc Woli 
1980h : 



( E V E x) 



(E x E* y ) 

(EyE* y ) 



I 
I 



Q 
u 



u 

-Q 



1 Italics are used to distinguish between the similarly notated Q band and Q Stokes parameter. 
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1 / 1 \ P / cos2 7 sin2 7 \ 

2 \ iy + 2Vsin27 -cos2 7 y (o) 

where we have set Stokes V = 0. The polarized component of the coherency matrix is a 
spin-two field on a sphere; the total power is the trace of the coherency matrix. 

The Crab Nebula [Tau A, 3C144, RA = 05 h 34 m 31 s , Dec=22°01' (J2000)] is the brightest 
polarized point source in the sky and provides a useful end-to-end check of the sign conven- 
tions and mapmaking pipeline. Figure [2] shows our measurement of the Crab in Q band 
(41 GHz) in /, Q, U, P, and 7. Note that its polarization direction (U ~ 0, Q negative), is 
perpendicular to the polarization of the Galaxy (U ~ 0, Q positive). The WMAP polariza- 
tion direction and intensity are in general agreement with previous measurements. Table [I] 
summarizes the results in all five frequency bands and previous measurements in our fre- 
quency range. A second check is needed to fully resolve the sign convention because with 
U = 0, 7 = ^pa- In Figure [2] we show that the polarization direction of the Centaurus A 
galaxy [Cen A, NGC5128 RA=1 3 h 25 m 27 s , Dec= -43°01'09" (J2000)] is consistent with that 



measured by iJunkes et al.l (119931 ) . 



Figures [3] and H] show the P and 7 maps of the full sky for all five frequency bands in 
Galactic coordinates. Figure shows a Lambert equal area projection of the Galactic polar 
region in K band. A number of features are immediately apparent to the eye. K band is 
strongly polarized over a large fraction of the sky, including the polar region. The North 
Polar Spur and its southern extension are clearly evident. The polarization has a coherent 
structure over large swaths of sky which translates into significant emission at low I. The 
polarization intensity decreases with increasing frequency but follows the same pattern. K 
band is a good monitor of polarized foreground emission as discussed below. Though not 
immediately apparent to the eye, there is somewhat more polarized emission at W band 
than V band. The uneven weighting due to the scan strategy is also evident as increased 



noise in the ecliptic plane (IBennett et al.ll2003bl . Figure 4). Figure M shows the K and Ka 
bands in Stokes Q and U. 

While foreground emission is visible with a high signal to noise ratio, the CMB polar- 
ization anisotropy is not, a situation unlike that for the temperature anisotropy. 



3. Systematic Errors 

Detection of the CMB polarization requires tight control of systematic errors, as small 
couplings to the temperature field or instrument will dominate the polarization signal. 
WMAP's differential nature and interlocked scan strategy suppress potential polarization 
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systematics in ways similar to the suppression for temperature systematics. The details are 
different however, and more complex because of the tensorial nature of the polarization field 
and the double difference required to measure the polarization. Throughout our analyses, 
th e overall level of syst em atic contaminat i on is assessed with null tests as described here and 



m 



Jarosik et all (120061 ) & iHinshaw et all (120061 ) . 



The mapmaking procedure is described in I Jarosik et al.l (120061 ). End-to-end simulations 
of the instrument and scan strategy, incorporating realistic models of the frequency response, 
foreground emission, and detector noise characteristics, are used to assess the possible lev- 
els of contamination. Interactions between the slow < 1 % drifts in the gain, non-uniform 
weighting across the sky, the 0.2% correlation due to the oppositely directed beams, the 
time series masking of the planets, and the 1/f noise are accounted for in the map solu- 
tion. In the following we discuss how the instrumental offset, gain/ calibration uncertainty, 
passband mismatch, main beam mismatch, polarization isolation and cross polarization, loss 
imbalance, and sidelobes affect the polarization maps. 

Offset and baseline drift — - The instrumental offset is the output of the detector in 
the absence of celestial signal. The average polarization offset in the Q, V, and W bands 
is 250 niK. Changes in this offset on time scales of minutes to hours arise from spacecraft 
temperature changes and from 1/f drifts in the amplifier gain acting on the 250 mK. To 
measure polarization at the level of 0.1 fiK, we require that changes in the baseline be 
suppressed by roughly a factor of 10 6 . The first step in achieving this is maintaining a stable 
instrument and environment. The physica l temperature of th e DAs averaged over a spin 
period changes by less than 5 parts in 10 6 (I Jarosik et al.l 120061 ) . suppressing changes in the 
baseline by a similar factor. T he second step in achieying this is through the baseline removal 
in the mapmaking algorithm (IHinshaw et al.ll2003bl ; I Jarosik et al.ll2006l ). 



If the precession of the satellite were stopped, the temperature data for £ > 1 would 
repeat in the time stream at the spin period (2.16 m). The offset, though, would change sign 
relative to the celestial signal at half the spin period enabling the differentiation of celestial 
and instrumental signals. (Alternatively, one may imagine observing a planet in which case 
the temperature data would change sign at half the spin period and an offset would be con- 
stant.) By contrast, with our choice of polarization orientations, the polarization data ATp, 
would repeat at half the spin period for some orientations of the satellite. Consequently, an 
instrumental offset would not change sign relative to a celestial signal upon a 180° space- 
craft rotation. Thus the polarization data are more sensitive to instrumental offsets than are 
the temperature data. In general, the polarization data enters the time stream in a more 
complex manner than does the temperature data. 



Calibration — An incorrect calibration between channels leads to a leakage of the tern- 
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perature signal into ATp, contam inating the p o larizat ion map. Calibration drifts cause a 
leakage that varies across the sky. iJarosik et al.l (l2003al ) show that calibration drifts on « 1 
day time scales are the result of sub-Kelvin changes in the amplifier's physical temperature. 
The calibration can be faithfully modeled by fitting to the physical temperature of each DA 
with a three parameter model. Here again WMAP's stability plays a key role. The residual 
calibration errors are at the m 0.2% level. These errors do not limit the polarization maps be- 



cause the bright Galacti c plane is mas 



Galactic latitude maps (IJarosik et al. 



still the first-year value, 0.5% (IJarosik et al.l 120061 ). 



red in the time ordered data when producing the high 
2006h. The overall absolute calibration uncertainty is 



Passband mismatch — The effective central frequencies (IJarosik et al.ll2003bl ; iPage et al. 



2003bl ) for AT X and AT 2 are not the same. This affects both the beam patterns, treated 
below, and the detected flux from a celestial source, treated in the following. The passbands 
for the A and B sides of one polarization channel in a DA may be treated as the same 
because the dominant contributions to the passband definition, the amplifiers and band 
defining filters, are common to both sides. 

Since WMAP is calibrated on the CMB dipole, the presence of a passband mismatch 
means that the response to radiation with a non-thermal spectrum is different fro m the 
response to radiation with a CMB spectrum (IKogut et al.ll2003l ; iHinshaw et al.ll2003bl ). This 
would be true even if the sky were unpolarized, the polarization offset zero, and the beams 
identical. The effect produces a response in the polarization data of the form: 



AT P = Ah - AI 2 + 

Q(n A ) cos 2-f A + U (h A ) sin 27,4 
-Q{n B ) cos 27b -U(n B ) sin 27^. 



(6) 



where AIi is the unpolarized temperature difference observed in radiometer one, and sim- 
ilarly for A/2- If these differ due to passband differences, the polarization data will have 
an output component that is independent of parallactic angle. Given sufficient paralactic 
coverage, such a term can be separated from Stokes Q and U in the mapmaking process. 
We model the polarized signal as Q cos 2j + U sin 27 + S* where the constant, S, absorbs the 
signal due to passba nd mismatch. We s olve for the mismatch term simultaneously with Q 
and U as outlined in IJarosik et al.l (120061 ). Note that we do not need to know the magnitude 
of the passband mismatch, it is fit for in the mapmaking process. The S map resembles 
a temperature map of the Galaxy but at a reduced amplitude of 3.5% in K band, 2.5% in 
the VI band, and on average ~ 1% for the other bands. The maps of S agree with the 
expectations based on the measured passband mismatch. 



Beamwidth mismatch — The beamwidths of each polarization on each of the A and B 
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sides are different. The difference between the A and B side beam shapes is due to the 
difference in shapes of the primary mirrors and is self consistently treated in the window 
function JPage et al.ll2003bh . The difference in beam shapes between A7\ and AT 2 is due to 



the mismatch in central frequencies^] 

This effect is most easily seen in the K-band observations of Jupiter. We denote the 
brightness temperature and solid angle of Jupiter with Tj and Oj, and the measured quanti- 
ties as Tj and Clj. Although the product TjVtj = TjClj is the same for the two polarizations 
(because Jupiter is almost a thermal source in K band), the beam solid angles differ by 



8.1% on the A-side and 6.5% on the B-side (IPage et al.ll2003al ). The primary effect of the 



beamwidth mismatch is to complicate the determination of the intrinsic polarization of point 
sources. 

The difference in beams also leads to a small difference in window functions between AT X 
and AT 2 . The signature would be leakage of power from the temperature anisotropy into the 
polarization signal at high t. We have analyzed the data for evidence of this effect and found 
it to be negligible. Additionally, as most of the CMB and foreground polarization signal 
comes from angular scales much larger than the beam, the difference in window functions 
can safely be ignored in this data set. 

Polarization isolation and cross polarization — Polarization isolation, X cp , and cross 
polarization are measures of the leakage of electric field from one polarization into the mea- 
surement of the orthogonal polarization. For example, if a source were fully polarized in the 
vertical direction with intensity I v and was measured to have intensity Ih = 0.0 II v with a 
horizontally polarized detector, one would say that the cross polar response (or isolation) 
is \X cp \ 2 = 1% or —20 dB. The term "polarization isolation" is usually applied to devices 
whereas "cross polarization" is applied to the optical response of the telescope. We treat 
these together as a cross-polar response. For WMAP, the off-axis design and imperfections 
in the orthomode transducers (OMT) lead to a small cross-polar response. The ratio of the 
maximum of the modeled crosspolar beam to the maximum of the modeled copolar beam is 
—25, —27, —30, —30, & —35 dB in K through W bands respectively. The determination of 
the feed and OMT polarization isolation is limited by component measurement. The maxi- 
mum values we find are: |A rp | 2 = —40, —30, —30, —27, & —25 dB for K through W bands 



respectively (IPage et al.ll2003bl ). We consider the combination of beams plus components 
below. 

Because WMAP measures only the difference in power from two polarizations, it mea- 



2 If the passbands were the same, the beam solid angles for ATi and AT 2 would be the same to < 0.5% 
accuracy. 
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sures only Stokes Q in a reference frame fixed to the radiometers, Quad- The sensitivity to 
celestial Stokes Q and U comes through multiple observations of a single pixel with different 
orientations of the satellite. The formalism that describes how cross polarization interacts 
with the observations is given in Appendix |A] To leading order, the effect of a simple 
cross polarization of the form X cp = Xe lY is to rotate some of the radiometer U into a Q 
component. The measured quantity becomes: 



AT P = Q^ d + Q B ad + 2X cos(Y)(U^ d + < d 



(7) 



where Q^ arf and Q B nd are the Stokes Q components for the A and B sides in the radiometer 
frame, similarly with U^ ad and t/J? ad . Note that in the frame of the radiometers Q Rad (Stokes 
Q in the B-side coordinate system) is — Qr^. This leads to the difference in sign conventions 
between the above and Equation [5j System measurements limit the magnitude of \X cp \ 2 
but do not directly give the phase, Y. Laboratory measurements of selected OMTs show 
Y = 90° ± 5°, indicating the effective cross polar contamination is negligible. 

We limit the net effect of the reflectors and OMT with measurements in the GEMAC 



antenna range ( IPage et al.ll2003bl ). We find that for a linearly polarized input, the ratios of 
the maximum to minimum responses of the OMTs are 1) —25, —27, —25, —25, —22 dB for K 
through W band respectively; 2) 90°±2° apart; and 3) within ±1.5° of the design orientation. 
Thus, we can limit any rotation of one compon ent into another to < 2°. The comparison 
of 7 derived from Tau A to the measurement by iFlett fc Henderson! (119791 ) in Table [1] gives 
further evidence that any possible rotation of the Stokes components is minimal. Based 
on these multiple checks, we treat the effects of optical cross polarization and incomplete 
polarization isolation as negligible. 

Loss imbalance — A certain amount of celestial radiation is lost to absorption by the 
optics and waveguide components. If the losses were equal for each of the four radiometer 
inputs their effect would be indistinguishable from a change in the gain calibration. However, 
small differ ences exist that prod uce a residual common-mode signal that is separable from the 
gain drifts (IJarosik et al.ll2003al ). The mean l oss difference (x\ m ) between the A- and B-s ides 
is accounted for in the mapmaking algorithm (IHinshaw et al.ll2003al ; IJarosik et al.ll2003al ) . In 
addition, the imbalance betw een the two polariza tions on a single side, the "loss imbalance 
imbalance," is also included (IJarosik et al.l 120061 ) . It contributes a term 2{L A T A + L B T B ) 
to ATp. Here T A,B is the sky temperature observed by the A,B side, and L A ' B is the loss 
imbalance betwe en the two polarizatio ns on the A, B side (see Appendix |Aj). The magnitude 
of L A ' B is < 1% J.Iarosik et al1l2003bh . 



A change in the loss across the bandpass due to, for example, the feed horns is a poten- 
tial systematic error that we do not quantify with the radiometer passband measurements 
(IJarosik et al.l l2003bl ). The magnitude of the effect is second order to the loss imbalance 
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which is 1%. We do not have a measurement of the effect. Nevertheless, as the effect mimics 
a passband mismatch, it is accounted for in the map solution. 

Sidelobes— Wh en the sidelobes corresponding to ATp are measured, there are two terms 
faarnes et alJl2003h . The largest term is due to the passband mismatch and is consistently 
treated in the mapmaking process. The second smaller term is due to the intrinsic polariza- 
tion. We assess the contribution of both terms by simulating the effects of scan pattern of the 
sidelobes on the Q and U polarization maps. The results are reported in lBarnes et al.l (120031 ) 
for the first-year polariza tion maps. In K band , the net rms contamination is 1 fiK outside 
of the KpO mask region (jBennett et al. 2003b ). The intrinsic polarized sidelobe pickup is 
< 1 /iK and is not accounted for in this three-year data release. The contamination is more 
than an order of magnitude smaller in the other bands. 



4. The Foreground Emission Model 



The microwave sky is polarized at all frequencies measured by WMAP. In K band 
the polarized flux exceeds the level of CMB polarization over the full sky. By contrast, 
unpolarized foreground emission dominates over the CMB only over w 20% of the sky. Near 
60 GHz and £ ~ 5, the foreground emission temperature is roughly a factor of two larger than 
the CMB polarization signal. Thus, the foreground emission must be subtracted before a 
cosmological analysis is done. While it is possible to make significant progress working with 
angular power spectra, we find that due to the correlations between foreground components, 
a pixel space subtraction is required. Table [2] gives the foreground emission levels in a region 
around the Galactic center. 

The two dominant components of diffuse polarized foregro und emission in the 23 — 



94 GH z range are synchrotron emission and thermal dust emission (jWeisslll984l ; iBennett et al. 
2003b ). Free-free emission is un polarized^l and spinning dust grains are expected to have 
polarization fractions of 1-2% (ILazarian fe Draind 120001) . The signal from polarized ra- 
dio sources is negligible (Table 9, iHinshaw et al.l 120061 ). The detected polarized sources 



are all well known, and among the brightest objects in the temperature source catalog. 
They include 3C273, 3C274 (M87, Vir A), 3C279, Fornax A, Pictor A, [HB93]2255-282, 
and [HB93] 0637-752 and are masked as discussed below. The potential impact of polarized 

man y 
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Synchrotron emission is produced by cosmic-ray electrons orbiting in the ~ 6 /iG total 
Galactic magnetic field. The unpolariz ed synchrotron comp onent has been well measured 
by WMAP in the 23 to 94 GHz range (IBennett et al.ll2003al ). The brightness temperature 
of the radiation is characteri zed by T(v) oc v^ a where the index —3.1 < (3 S < —2.5 varies 



considerably across the sky (IReich fc Reich! Il988l ; lLawson et al.l 119871 ) . In the microwave 
range, the spectrum reddens {(3 a tends to more negative values) as the frequency increases 
faandav fc Wolfenddelll99ll ). 



Synchrotron radiati on can be strongly polarize d in the direction perpendicular to the 
Galactic magnetic field (jRybicki fc Lightmanl Il979l ). The polarization has been measured 
at a n umber of frequencies [from Leiden between 408 M Hz to 1.4 GHz (IBrouw fc Spoelstra 



19761 : IWolleben et al.l 120051 ). from Parkes at 2.4 GH z flDuncan et al 



199 



5. 



1999|), and by 



the Medium Galactic Latitude Survey at 1.4 GHz (jUyamker et al.l Il999l )]. At these low 
frequencies, Faraday rotation alters the polarization. Electrons in the Galactic magnetic field 
rotate the plane of polarization because the constituent left and right circular polarizations 
propagate with different velocities in the medium. In the interstellar medium, the rotation 
is a function of electron density, n e , and the component of the Galactic magnetic field along 
the line of sight, B\\, 



A9 = 420 c 



1 GHz\ 2 f L/1 kpc 



v 



dr 



n. 



0.1 CHI" 



B 



1 fiGJ 



where the integral is over the line of sight. With n e ~ 0.1 cm" 3 , L ~ 1 kpc, and Bu ~ 1 fiG, 
the net rotation is A8 ~ 420°/z/ 2 , with v in GHz. At WMAP frequencies the rotation is 
negligible, though the extrapolation of low frequency polarization measurements to WMAP 
frequencies can be problematic. In addition there may be both o bservational and astrophys- 
ical depolarization effects that are different at lower frequencies ( iBurnl Il966l : ICioffi &: Jones 
19801 ; ICortiglioni fc Spoelstralll995l ). Thus, our method for subtracting the foreground emis- 
sion is based, to the extent possible, on the polarization directions measured by WMAP 



The other dominant component of polarized foreground emission comes from thermal 
dust. Nonspherical dust grains align their long axe s perpendicularly to the Galactic mag- 
netic field through the Davis-Greenstein mechanism (IDavis fc Greensteinlll95ll ). The aligned 
grains preferentially absorb the component of starlight polarized along their longest axis. 
Thus, when we observe starlight we see it polarized in the same direction as the magnetic 
field. These same grains emit thermal radiation preferentially polarized along their longest 
axis, perpendicular to the Galactic magnetic field. Thus we expect to observe thermal dust 
emission and synchrotron emission polarized in the same direction, while starlight is polarized 
perpendicularly to both. 
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In Section §4.1[ we describe a physical model of the polarized microwave emission from 
our Galaxy that explains the general features of the WMAP polarization maps. However 
this model is not directly used to define the polarization mask or to clean the polarization 
maps. We go on to define the polarization masks in §4.21 and in §4.31 we describe how we 
subtract the polarized foreground emission. 



4.1. The Galaxy Magnetic Field and a Model of Foreground Emission. 



In the following, we present a general model of polarized foreground emission based on 
WMAP observations. We view this as a starting point aimed at understanding the gross 
features of the WMAP data. A more detailed model that includes the wide variety of external 
data sets that relate to polarization is beyond the scope of this paper. 

For both synchrotron and dust emission, the Galactic magnetic field breaks the spatial 
isotropy thereby leading to polarization. Thus, to physically model the polarized foreground 
emission we need a model of the Galactic magnetic field. As a first step, we note that the 
K-band polarization maps suggest a large coherence scale for the Galactic magnetic field, as 
shown in Figure [31 

We can fit the large-scale field structure seen in the K-band maps with a gas of cosmic 
ray electrons interacting with a magne tic field that follows the sp i ral arms. The Galactic 
magn etic field can be quite complicated (|Beck!l2006l ; lHan et al.ll2006l ; lReichll2006l ; IWielebinski 
20051 ): there are field direction reversals in the Galactic pla ne; the field strength dep ends 
on length scale, appearing turbulent on scales < 80 kpc (e.g.. iMitner fc Spanglerlll996l ) ; and 
the fi eld strength of the large-scale field depends on the Galactocentric radius (e.g., [Beck 
20011 ) . Nevertheles s, most external spiral galaxie s have magnetic fields th at follow the spiral 
arm pattern (e.g., IWielebinskil 120051 ; iBeck et al.lll996l ; ISofue et al.lll986l ). Inspired by this, 
we model the field in cylindrical coordinates as: 



B(r, 0, z) — B Q [ cos ip(r) cos x{z)r + 
sin ip(r) cosx{z)4> + 
sin x(z)z] 



(9) 



where ip(r) = ipo + ^iln(r/8 kpc), x{ z ) = Xotanh(z/l kpc), r and z are measured in kpc 
with respect to the center of the Galaxy, r ra nges from 3 kpc to 20 kpc, and the angles 
are in degrees. The coordinates follow those in lTaylor &: Cordesl ( 119931 ) . For a fixed radius, 
— B — has the same value at all azimuths. We term the expression the logarithmic spiral 
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arm (LSA) model to disti nguish it from previou s forms. We take 8 kpc a s the distance to 
the center of the Galaxy (jEisenhauer et al.ll2003l ; iReid &: Brunthalerll2005l ). The values are 
determined by fitting to the K-band field directions. While the tilt, x( z ) with xo — 25°, and 
the radial dependence, tp(r) with ipi = 0?9, optimize the fit, the key parameter is ipo, the 
opening angle of the spiral arms. We find that the magnetic field is a loosely wound spiral 
with ?/>o — 35°. 

To model the cosmic r ay electrons, we assume the y have a power-law distribution with 
slope @p = — (2f3 s + 3) = 3 ( Rybicki fc Lightman 1979 ) and are distributed in a exponential 
disk with a scale height of ha = 1 kpc and a radial scale length of h r = 5 kpc (e.g., 
Drimmel &: Spergelll200ll ) as 



n. 



no exp(— r / h r )sech 2 (z / h 



(10) 



While the amplitude of the signal is sensitive to the details of the cosmic ray distribution 
and the magnetic field structure, we may estimate its overall structure with the smooth field 
model (Eq. [U]) and cosmic ray distribution. The model predictions are not very sensitive to 
the assumed scale height and scale length. We compute the polarization direction in this 
simple model as: 



n U(n) 
tan2 T n = 

Q[n) 



J n e (x,h)2B s (x,n)B t (x,n) dx 



f n e (x, n) [B^(x, n) — B\[x, n)] dx 

where h is the line-of-sight direction, x is the distance along that direction, n e is the elec- 
tron distribution described above, and B t and B s are orthogonal components of the field 
perpendicular to the line of sight, with B t the component perpendicular to the z axis of the 
Galactic plane. The parameters of the LSA model are determined by fitting the predicted 
directions, Equation [TTJ to the measured the K-band field directions. 

Figure [7] shows the predicted magnetic field orientation for the LSA model. The actual 
direction has a 180° ambiguity. In the plane, the field lines are parallel to the Galactic 
plane and the polarization projects into positive Stokes Q. Near the Galactic pole, the field 
lines point along the spiral arm direction. When projected into Q and U, this leads to 
7 rotating around the pole. We assess the agreement between the model field orientation 
and the orientation inferred from the K-band polarization with the correlation coefficient 



iBennett et al.l (|2003a[) uses 7 in place of p 
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r = cos(2(7 mo< iez — Jdata)), and take the rms average over 74.3% of the sky (outside the P06 
mask described below). For our simple model the agreement is clear: r = 0.76 for K band. 

Using using rotation measures derived from pulsar observations in t he Ga l actic p lane, 
one fi nds instead a spiral arm opening angle of ipo ~ 8° as reviewed in iBeck! ( 120011 ); lHan 
( 120061 ). Our method is more sensitive to the fields above and below the plane; and, unlike 
the case with pulsars, we have no depth information. It has been suggested by R. Beck 
and others that the north polar spur may drive our best fit value to ipo — 35°. Though 
the agreement between our simple model and the K-band polarization directions indicate 
that we understand the basic mechanism, more modeling is needed to connect the WMAP 
observations to other measures of the magnetic field. 

For a power law distribution of electrons moving i n a homogeneous magnetic field, the 
polarization fraction is Il s = (p + l)/(p + 7/3) ~ 0.75 (IRybicki &: Lightmanlll979l ). Because 
the field direction changes as one integrates along the line of sight, there is a geometric sup- 
pression of the amplitude of the polarization signal. We estimate this geometric suppression 

as 

P(n) 

9syne(n) = ^ ~TT^ ; (12) 



where all quantities are determined from the model: P(n) = a/ Q 2 + U 2 and I is found by 
integrating the perpendicular component of the magnetic field, (B 2 + B 2 ) 1 ^ 2 , and cosmic ray 
distribution al ong the line o f sight . The result is shown in Figure Similar results have 
been found by lEnBlin et al.l (120061 ). This geometric reduction factor ranges from unity to 



zero. 



4- 1.1. Comparison to Low Frequency Observations 



The polarization o f a nu mber of spiral galaxies similar to the Milky Way has been mea- 
sured by lDumke et al.l (If 9951 ). The observations are at 10.55 GHz and thus probe primarily 
synchrotron emission. For one of the best measured edge-on spirals, NGC 891, they find: (1) 
at distances ~ 5 kpc off the galactic plane the polarization fraction can be >10%; and (2) 
in the pla ne, at heights < 0.5 kpc, the polarization fraction drops to <5%. Simil ar behavior 
is seen by ISukumar fc Allenl (119911 ) at 5 GHz. In addition, iHummel et al.l (1199 ll ) show that 
(3) between 0.66 GHz and 1.5 GHz the spectral index ranges from f3 s = —2.5 in the plane 
to j3 s = —3.5 well off the plane. WMAP observes qualitatively similar behavior in K band. 



At 408 MHz, lHaslam et al.l (Il982l ) have surveyed the Galactic plane in intensity. At 
this frequency, synchrotron emission dominates maps. We test the magnetic field model by 
extrapolating the 408 MHz measurements to 22 GHz (an extrapolation of 40 in frequency 



-17- 



and over 10,000 in amplitude): 

Qmodel{h) = qluasin) 



Umodelin) 



qlHas{n) 



22 
0.408 

22 
0.408 



^s9sync{n) cos(27 mode j) 



Ti s gsync{n) sm(2^ mode i) 



(13) 



where q is the ratio of the homogeneous field strength to the total field strength. Note that 
the model effectively has only one free parameter: an overall amplitude, which is described 
by a degenerate combination of the spectral index, (3 S and q. For (3 S = —2.7, the best fit 
value for q is 0.7. This implies that the energy in the large scale field is roughly the same as 
the energy in sm all scale fields, consistent with other results for the Milky Way (jjones et al. 
19921 : iBeckl l200ll . and references therein) . 



Figure M compares the K-band polarization signal to the extrapolated 408 MHz maps. 
Given the simplicity of the model (uniform cosmic ray spectral index, p, and a uniform 
LSA field), the agreement i s remarkably good. T he largest deviations are seen near spiral 
arms. Recent observations (jEnomoto et all 120021 ) suggest that cosmic rays are accelerated 
in star-forming regions. If most cosmic rays are accelerated in spiral arms and then diffuse 
away from the arms, we would expect a flatter spectral index in the arms, consistent with 
the observations. In Figure [TO] we show that the radio loops (jBerkhuijsen et all Il97ll ) seen 
at 408 MHz, probably from supernovae or "blowouts," are also seen in the WMAP data. 



4-1-2. Starlight Polarization and Polarized Dust Emission 



Measurements of starl ight polarization ser v e as a template for th e analysis of polarized 



microwave dust emission flFosalba et al.1 120021 : iBernardi et a 



several catalogs of optical dust po l arization measurement s "dHeilesI l2000l : iBerdyugin et al.l 



2003]). We have combined 



200ll : IBerdyugin fc Teerikorpil 120021 : IBerdyugin et al.l 120041 ) to construct a template for the 
magnetic field direction in dusty environments. Since there are significant variations in the 
dust column density, we only use the mea sured direction to construct the dust template. The 
dust layer has a scale height of 100 pc (IBerdyugin fc Teerikorpil 1200 ll : iDrimmel fc Spergel 
200ll ) . Observations toward the Galactic poles suggest that most of the dust absorption 
occurs within 200 pc. To select stars outside the dust column for |6| > 10°, we limit the 
sample to the 1578 stars with heliocentric distances greater than 500 pc. For |6| < 10°, the 
model is problematic because there is ample dust emission from distances further away than 
the stars sample. 



We represent the starlight polarization data, (Q*, U+), in terms of a polarization ampli- 
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tude, P+ and direction, 7*: 

Q+ = cos(27^) 

U+ = P*sin(2 7 *) (14) 

We then smooth the starlight data by convolving (Q+/P+) and (U+/P+) with a Gaussian 
window with a FWHM of 9?2. The smoothing is required because the measurements are 
coarsely distributed. As a result, this dust model is applicable only for £ < 15 and \b\ > 10°. 
Above, 7* describes the direction of this smoothed starlight polarization field. We can 
quantify the agreement between the starlight and WMAP K-band polarization measurements 
by computing their correlation in each pixel, Z = cos(2(7*— 7^)+7r) where jk is the direction 
in K band. Figure [11] shows a plot of the correlation as a function of position. The median 
correlation coefficient is 0.72 implying that the dust and K-band directions typically agree 
to 20°. Because of noise in both the K-band and starlight maps, this is an underestimate 
of the correlation. Nevertheless, the correlation tells us that the basic model relating the 
starlight, the dust, synchrotron emission, and the magnetic field agrees with observations. 



4- 1.3. Thermal Dust Emission 



Based on the detection of starlight polarization, thermal dust emission is expected to be 
polarized at millimeter and s ub-millimeter waye lengths. Archaeops has detected polarized 
thermal emission at 353 GHz (IBenoit et al.ll2004l ). An extrapolation from this high frequency 
suggests that WMAP should see polarized thermal dust emission at 94 GHz. Here, we report 
on the WMAP detection of dust polarization at 94 GHz. 

We generate a template for the dust polarization by using the Maximum Entropy 
Method (MEMjB dust intensity map (IBennett et al.ll2003al ). the smoothed polarization di- 
rection from the starlight, and the model geometric factor for the dust layer: 



Qdust{v) — I ( 
U d ust{v) = 



r dust(n)Udgdu 



lustifi) cos{2<~i dust 
I d ust(n)U d g dust (n) sin(27 



dust J 



(15) 



where jdust — 7* + 7r /2 is the smoothed starlight polarization direction. The geometric 
suppression factor for the dust, gdust, is computed along the same lines as g sy nc in Equation 



5 The dust, free-free, synchrotron, and CMB MEM maps are derived from a maximum entropy solution 
to the five WMA P bands, the FDS dust map (jFinkbeiner et al.lll999| ), the Haslam map, and a Ha map 
|Finkbeinerll2003h 
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and is shown in Figure To compute 1(h) for the analog to Equation [T2J we assume the 
dust has a scale height of 100 pc and a radial scale length of 3 kpc. To find P(h) we use the 
LSA magnetic field model. The fractional polarization, 11^ = 0.05, is found with a best fit of 



the m odel to the data. Similar results are found using the FDS dust map (IFinkbeiner et al. 



19991 ) instead of the MEM dust map. The uncertainty is estimated to be 50%. 



Figure [12] compares this predicted pattern of polarization to the cleaned W-band ob- 
servations. We use the K-band synchrotron template to clean Q, V and W bands and then 
use the Q and V band maps to remove the CMB polarization signal from the W-band maps; 
though removing the CMB component is not necessary. The W-band map is then smoothed 
with a 10° beam for plotting. T he appearance of the dust polarization signal pattern is sim- 
ilar to that found by Archeops (IPonthieu et al.ll2005l . Figures 2 & 3). However, the signal to 
noise ratio is low due to the low level of polarized dust emission at 94 GHz. The predominant 
feature is that the plane is dominated by positive Stokes Q emission. A visual comparison to 
the model is less robust. One must keep in mind that since stars are heavily obscured in the 
plane, the model is not expected to be accurate in the plane. Nevertheless, since Stokes Q 
emission corresponds to the dominant horizontal magnetic field, one does not have to sample 
too deeply to pick it up. Similarly, we interpret the poor correlation between the model U 
and the observed U as due to the insufficient sampling of other magnetic field directions by 
rather limited depth of the stars. Some common features between the model and W-band 
data are seen for \b\ > 10°. Fits of the data to the model are given in Section §4.31 Clearly, 
more integration time and more stellar polarization measurements are needed to fill out the 
model. 



4- 1-4- Spinning Dust Emission? 



Electric dipole emission from rapidl y spinning dust grains is potentially a significant 
source of emission at WMAP frequencies ( IEricksonlll957l ; lDraine &: Lazarianlll998l ). Thermal 
fluctuations in the magnetization of magnetize d grains may also be a p o tentially significant 
sourc e of emission at microwave wavelengths (jDraine fc Lazarianl 1 19991 ; iPrunet fc Lazarian 
19991 ). Both have been proposed as an explanation for the correlations seen between thermal 
dust em ission at 140 /xm a nd micro wave emission in m any cosmic b ackground experiments: 



COB E (jKogut et alJll996h . OVRO flLeitch et all 1 1 99 7Q Saskatoon f de Oliveira-Costa et al. 



1997h the 19 GHz Sur vey ddeOliveira-Costa et al.lll998h . Tenerife ( 



de Oliveira-Costa et al. 



19991 l20Q4h . Pvthon V jMukheriee et al.l 2003h. and COSMOSOMAS iFernandez-Cerezo et al 
2006J). 



The spectral shape of spinning dust emission can be similar to synchrotron emission 
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i n the 20-40 GHz ran ge. Thus models with either variable synchrotron spectral index 
(IBennett et al.l l2003br) or with a spinni ng dust spectrum with a suitably fit cutoff fre- 
quency ( Lagache 12003 ; Finkbeiner 2004 ) can give reasonable fits t o the data. However, 



at v < 20 GHz there is a considerable body of evidence, reviewed in IBennett et al.l (l2003bl ) 
and iHinshaw et al.l (120061 ). that shows (1) that the synchrotron index varies across the sky 
steepening with increasing galactic latitude (as is also seen in WMAP ) and (2) that in other 
galaxies and our galaxy there is a strong correlation between 5 GHz synchrotron emission 
and 100 jum (3000 GHz) dust emission. The combination of these two observations imply 



that the v < 40 GHz WMAP 



oregr ound emission is dominated by synchrotron emission 



as discussed in IHinshaw et al.l (120061 ). Nevertheless, we must consider spinning dust as a 



possible emission source. While on a Galactic scale it appears to be sub-dominant, it may 
be dominant or a significant fraction of the emission in some regions or clouds. 

Spinning dust models predict an unambiguous signature in intensity maps: a t 5-15 GHz, 
the du st e mission should be signifi c antly less than the synchrotron emission. iFinkbeiner 
(120041 ) and Ide Oliveira-Costa et al.l (120041 ) argue that the Tenerife and Green Bank data 
show evidence for a rising spectrum between 10 and 15 GHz, suggesting the presence of 
spinning dust. Observations of indiv i dual compact clouds also show evidence for spinning 
dust emission (IFinkbeiner et al.ll2002l . 120041 ; IWatson et al.ll2005l ) t hough the signature i s not 
ubiquitous. The status of the observations is discussed further in IHinshaw et al.l (120061 ). 



The WMAP polarization measurements potentially give us a new way to distinguish 
between synchrotron and dust emission at microwave frequencies. While synchrotron emis- 
sion is expected to be highly polarized, emission from spinning dust grains is thought to 
be weakly polarized. While promising, the signature is no t unique as a tangle of m agnetic 
field lines can also lead to a low polarization component (jSukumar fc Allenl Il99ll ) as seen 
at 5 GHz where spinning dust emission is expected to be negligible. Using a model for the 
polarization fraction of the synchrotron emission based on the LSA structure, we separate 
the microwave intensity emission into a high and low polarization component: 

P v (h)/qU s g sync (n) 

I u (n) ~ - ISS B ~ ij?"'" (16) 

where I v and P v are the intensity and polarization maps at frequency v. For notational 
convenience, we use v = K, Ka, Q, V, W. Iqmb l § ^ ne foreground-free CMB map made with 
a linear combination of WMAP bands, and Ipp M ' u is the MEM free-free map for band 



v (IBennett et al.l l2003bl ). In effect, we use the WMAP polarization maps to extract the 
intensity map of the low-polarization component in the data. 

Figure [13] co mpares the morpholog y of the low polarization K-band map to the W-band 



MEM dust map (IBennett et al. 



2003al ). Even in this simple model based on a number of 
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assumptions, the agreement in morphology is striking. We quantify this by computing the 
rms deviation between the two scaled maps, 

d 2 = £(Mn) - QJ&W (17) 

where W is the W-band map, the scale factor is a = 0.105, and the sum is taken over pixels. 
We find d = 0.05. In other words, we can "predict" the distribution of dust in W band from 
just the K band intensity and polarization maps. The low polarization fraction component 
has a spectral index of /3 — —2.6 between K and Q bands. This correlation between the 
low polarization emission regions at 22-45 GHz and the thermal emission at 90 GHz and 
higher may be interpreted as either a very tight correlation between tangled field lines in 
star forming (dusty) regions or as evidence for spinning dust emission. More polarization 
data, v < 22 GHz observations, and extensive modeling are needed to conclusively delineate 
the magnitude and morphology of the various components. 



4.2. Masking Polarized Foreground Emission 

To compute the CMB power spectrum, we must mask the regions with the brightest 
foreground emission. For polarization we create a set of m asks with a process th at is some- 
what analogous to the creation of the temperature masks (IBennett et al.ll2003bl ). First, the 
K-band Q and U polarization maps are use d to compute a posi tive-definite HEALpix r4 § P 
map. From this a noise-bias variance map (IJarosik et al.ll2006l ) is subtracted. The rectified 
noise-bias correction is small because of the coarse resolution at r4. A histogram of pixel 
polarization amplitudes in this noise-bias-corrected map approximates a power law. The 
peak is near the zero pixel value, there are just a few negative pixels (due to the noise bias 
correction), and there is a long positive tail. 

Unlike the process in which the temperature masks were created, there is no natural 
cut level based on the histogram peak. Instead, the cuts are given in terms of the mean of 
the noise bias corrected map of P at K band. The cut level at the mean is denoted "P10". 
The cut level at 0.2 times the mean is "P02" , etc. For each cut level, a preliminary mask is 
made by setting r4 pixels greater than the cut level to 1, and all others to 0. This mask is 
expanded to r9 and smoothed by a 7?5 FWHM Gaussian. This mask map is set back to all 
0s and Is using the 0.5 level as a cut-off and the sense of the mask is reversed, so that the 



6 The n umber of pixels is 12N ? irl „ where N s u e 
notation in lBennett et al.l (|2003bl ). 



2 4 for r4, or resolution 4 ( Gorski et alJ ll998). See 
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masked-out parts of the sky have zeros (the WMAP convention). The above process results 
in a synchrotron polarization mask. 

In the case of temperature masks, we found that additional masking based on the higher 
frequency bands was redundant. This is not the case with polarization. Thus we make a 
dust polarization mask in a similar manner. We begin with the first-year MEM dust model 
box-averaged to r7. Half the maximum value found in a subset of pixels in the polar caps 

> 60°) is adopted as the cut-off level. A preliminary mask is made by setting r7 pixels 
greater than the cut-off level to unity, and all others to zero. This mask is then resolution 
expanded to a r9 map, smoothed by a 4?0 Gaussian, and set back to digital levels with a 
0.5 cut-off. The sense of the mask is then reversed to fit the WMAP convention. Each 
synchrotron polarization mask is ANDed with the (constant) dust polarization mask and a 
constant polarized source mask. 

We find, in general, that th e extragalactic point sources are minimally polarized in the 



WMAP bands, as discussed in iHinshaw et al.l (120061 ). We construct a source mask based 



on the exceptions. The most significant exception (not already covered by the synchrotron 
or dust polarization masks) is Centaurus A, an extended and polarized source. We found 
excellent agreement between WMAP and previously published maps of Cen A (Figure [2]). 
Based on this information, we custom-masked the full extent of Cen A. Six other bright 
polarized sources that we masked are Fornax A, Pictor A, 3C273, 3C274, 3C279, PKS 1209- 
52. (Some bright polarized sources already covered by the synchrotron and dust mask regions 
include: 3C58, Orion A, Taurus A, IC443, 1209-52, W51, W63, HB21, CTB104A). We have 
determined that, for most applications, the mask that we call "P06" is the best compromise 
between maximizing usable sky area while minimizing foreground contamination. With the 
above considerations, the P06 mask masks 25.7% of the sky, mostly near the Galactic plane. 
We use the terminology "outside the P06 mask" to refer to data in the 74.3% of the sky left 
for cosmological analysis. Various masks are shown in Figure [T4] 



4.3. Removing the Polarized Foreground Emission from the Maps 

Based on our analysis of the Galactic foreground emission, we have generated syn- 
chrotron and dust template maps for the purposes of foreground removal. The template 
maps are fit and subtracted from the Ka through W band data to generate cleaned maps 
that are used for CMB analysis. We assess the efficacy of the subtraction with x 2 an d by 
examining the residuals as a function of frequency and multipole £, as described in §5.21 



We use the K-band data to trace the synchrotron emission, taking care to account for the 
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(relatively weak) CMB signal in the K-band map when fitting and subtracting the template. 
For dust emission, we construct a template following Equation (fT5j) tha t is based on the 



starli ght- derived polarization directions and the FDS dust model eight (IFinkbeiner et al. 



19991 ) evaluated at 94 GHz to trace the dust intensity, /dust- We call this combination of 



foreground templates "KD3Pol" . 

The synchrotron and dust templates are fit simultaneously in Stokes Q and U to three- 
years maps in Ka through W bands. The three-year maps are constructed by optimally 
combining the single-year maps for each DA in a frequency band. Specifically 

[Q V ,U V ] = (^Nrr^Nr 1 ^,^] (18) 

i i 

where i is a combin ed year and DA index, [Qi, Uj\ is a polarization map degraded to r4 



( Jarosik et al.l 120061 ). and N f 1 is the inverse noise matrix for polarization map i. The fit 
coefficients, a s and are obtained by minimizing x 2 , defined as 

2 _ \- {[Qu, Uv] - a s AQs, Us] - adAQd, Ud]f hq s 

where [Q S ,U S ] is the K-band polarization map (the synchrotron template), [Qd,Ud] is the 
dust template, and [cr^cx 2 /] is the noise per pixel per Stokes parameter in the three-year 
combined maps. We have tried using optimal (N _1 ) weighting for the fits as well, and found 
similar results for the coefficients. The results reported here are based on the simpler diagona l 



weighting. The fit is evaluated for all pixels outside the processing mask (j Jarosik et al.ll2006l ). 



The fit coefficients are given in the top half of Table [3j For each emission component 
we also report the effective spectral index derived from the fit: (3 s {vk-iV) for synchrotron 
emission, and @d{v, vw) f° r dust. These results indicate that the spectrum of the component 
traced by K-band is systematically flattening with increasing frequency, which is unexpected 
for synchrotron emission. This behavior is statistically significant, and is robust to variations 
in the dust model and the data weighting. We do not have a definitive explanation for this 
behavior. 

To guard against the possibility of subtracting CMB signal, we modified the template 
model as follows. We take the four synchrotron coefficients in Table [3] and fit them to a 
spectrum model of the form 

<x,(v) = ot sfi - g{v){v/v K y° +a c , (20) 

where a Sj o, /3 S , and a c are model parameters that are fit to the a s (i/), and g(y) is the 
conversion from antenna temperature to thermodynamic temperature at frequency v. This 
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results in a modified set of synchrotron coefficients that are forced to follow a power-law that 
is largely determined by the Ka and Q-band results. Specifically, the modified coefficients 
are given in Table [3J The implied synchrotron spectrum is f3 s = —3.33. This results in a 
12% reduction in the synchrotron coefficients at Q-band, and a 33% reduction at V-band. 
However, because the K-band template is dominated by an i = 2 E-mode signal (see §5.ip . 
this change has a negligible effect on our cosmological conclusions, which are dominated by 
E-mode signal at i > 2. A comparison of selected "before and after" cleaning maps is shown 
in Figure [T51 

We also account for the cleaning in the map error bars. Since the K-band data are a 
combination of synchrotron and CMB emission, subtracting a scaled version of K band from 
a higher frequency channel also subtracts some CMB signal. If the fit coefficient to the higher 
frequency channel is a , then the cleaned map is M'iv) = (M(y) — a^Miv = K))/(l — ao), 
where M is the map and v denotes the frequency band. The maps we use for cosmological 
analysis were cleaned using the coefficients in the bottom half of Table [31 The factor of 
1/(1 — ao) dilates the noise in the cleaned map. To account for this in the error budget 
we scale the covariance matrix of the cleaned map by a factor of 1/(1 — ao) 2 . Additionally, 
we modify the form of the inverse covariance matrix by projecting from it any mode that 
has the K-band polarization pattern: N -1 — > N' -1 , where N' _1 [Qx, Uk] = 0. This ensures 
that any residual signal traced by K-band (due, for example, to errors in the form of the 
spectrum) will not contribute to cosmological parameter constraints. 

One measure of the efficacy of the foreground removal is the change in x 2 , relative to 
a null signal, between pre-cleaned and cleaned maps. Table H] gives the values for the full 
sky and the P06 cut. In both cases the full pixel covariance matrix was used to compute \ 2 
for Stokes Q and U simultaneously. For the full sky the number of degrees of freedom, u, is 
6144 (twice the number of pixels in an r4 map) and outside the P06 mask v = 4534. Note 
the large A% 2 achieved with just a two parameter fit. By comparing the full sky to the P06 
X 2 , we find that the starlight-based dust template is insufficient in the plane as discussed 
in §4.1.31 We also see that outside the P06 mask, that Q and V bands are the cleanest 
maps and that they are cleaned to similar levels. Since x 2 jv for Q and V bands is so close 
to unity for the cleaned maps, it is no longer an effective measure of cleaning. Instead, we 
examine the power spectra t by t to assess the cleaning, and then test the sensitivity of the 
cosmological conclusions to cleaning by including Ka and W band data. 

We have tried a number of variants on the KD3Pol cleaning. We find, for example, 
that setting g dust = 1 across the sky has negligible effect on the fits or the derived optical 
depth. Alternatively, when one uses the K-band polarization direction to trace the dust 
direction, 7d us t = 7k in Equation [T5| the cleaning is not as effective. Outside the P06 mask, 
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the reduced x 2 in the Q and V band maps is 1.022 and 1.016 as compared to 1.014 and 
1.018 for the starlight-based directions. Thus the latter are used. Regardless of template, we 
find that our cosmological conclusions are relatively insensitive to the details as indicated in 
Figure [26J 



5. Power Spectra 

The Q & U maps are well suited to analyzing foreground emission, are useful for com- 
paring to other polarization maps, and have straightforward noise properties. However, 
they are not well suited to quantifying the CMB polarization anisotropy because their 
definition is coordinate dependent. The Q and U maps may be transformed into scalar 



and pseudo-scalar quantiti es called E and B modes (ISeljaklll997l ; iKamionkowski et al.lll997l ; 
Zaldarriaga fc Seljald 119971 ) . E and B are so named because they comprise a curl- free and 
divergence-free decomposition of the spin-2 polarization field, analogous to static electric and 
magnetic fields. The problem of separating E and B mod es with an unevenly s ampled and 



cut sky has been consid ered by a number of authors (e.g., iTegmark et al.ll2000l ; iLewis et al. 
20021 ; iBunn et al.l 120031 ). In our analysis, we work directly with Q and U map s to produce 
the E and B angular power spectra. The conventions follow Appendix A of iKogut et al. 



Fundamental symmetries in the production and growth of the polarized signal select 
the possible configurations for the CMB polarization. Scalar (density) perturbations to the 
matter power spectrum give rise to T and E modes. Tensor perturbations (gravitational 
waves) give rise to T, E, and B modes primarily at £ < 20CQ Both scalar and tensor 
perturbations can produce polarization patterns in both the decoupling and reionization 
epochs. Vector perturbations (both inside and outside the horizon) are redshifted away 
with the expansion of the universe, unless there are active sources creating the vector modes, 
such as topological defects. We do not consider these modes here. 

At the noise levels achievable with WMAP , the standard cosmological model predicts 
that only the E mode of the CMB polarization and its correlation with T will be detected. 



7 In this paper we do not use the rotationally invariant Q' and U' of IKogut et al 



(200 



3). 



8 For r < 0.03 and I > 70, primordial B modes are dominated by the gravitational lensing of E modes. 

9 Vector modes are produced by purely rotational fluid flow. Based o n the fit of the adia batic ACDM 
model to WMAP TT data, the contribution of such modes is not large ( Spergel et al. 2003). However, a 
formal search for them has not been done. 
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The B-mode polarization signal is expected to be too weak for WMAP to detect, while the 
correlations of T and E with B is zero by parity. Thus the TB and EB signals serve as a 
useful null check for systematic effects. The polarization of foreground emission is produced 
by different mechanisms. Foreground emission can have any mixture of E and B modes, it 
can be circularly polarized (unlike the CMB), and E and B can be correlated with T. 

We quantify the CMB polarization anisotropy with the Cj E , Cf E , and Cf B angular 
power spectra, where 



C: 



XY 



i X Y*\ 
\ a lm a lml- 



(21) 



Here the "()" denote an ensemble average, aj m are the multipoles of the temperature map, 
and af m , af m are related to the spin-2 decomposition of the polarization maps 



[Q±iU](x) = J2J2 



(22) 
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(IZaldarriaga fc Seljaki 119971 ). The remaining polarization spectrum combinations (TB, EB) 
have no expected cosmological signal because of the statistical isotropy of the universe. 

We compute the angular power spectrum after applying the P06 polarization mask using 
two methods dependin g on the £ range. A ll power spectra are initially based on the single- 



year r9 Q and U maps ( Jarosik et al. 2006). For I > 23 B we compute the po wer spectrum 



following the method outlined in iHivon et al.l (120021 ). and lKogut et al.l (120031 . Appendix A) 
as updated in iHinshaw et al.l (120061 ) and Appen dix IB.2I The statistical weight per pixe l 



is N b s /<TQ where o"o is the noise per observation (j Jarosik et al.ll2006l ; IHinshaw et al.l 120061 ) . 
Here N b s is a 2x2 weight matrix that multiplies the vector [Q, U] in each pixel 



N, 



obs 



Nq N qu 
N QU N v 



(24) 



where Nq, Njj, and Nqu are the elements of the weight arrays provided with the sky map 
data. Note that the correlation between Q and U within each pixel is accounted for. We 
refer to this as "iV & s weighting." From these maps, only cross power spectra between DAs 



10/ 



= 23 = 3N S ide — 1 is the Nyquist limit on I. For some analysis methods ( CDj) we use HEALPix r3 for 
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and years are used. The cross spectra have the advantage that only signals common to two 
independent maps contribute and there are no noise biases to subtract as there are for the 
auto power spectra. The covariance matrices for the various Cg are given in Appendix IC.3I 

For I < 23 we mask and degrade the r9 maps to r4 (see the last paragraph of Appendix D 



and I Jarosik et al.l 120061 ) so that we may use the full r4 inverse pixel noise matrix, N _1 , to 
optimally weight the maps prior to evaluating the pseudo-CV This is necessary because the 
maps have correlated noise that is significant compared to the faint CMB signal. By "N _1 
weighting" the maps, we efficiently suppress modes in the sky that are poorly measured 
given the WMAP beam separation and scan strategy (mostly modes with structure in the 
ecliptic plane). We propagate the full noise errors through to the Fisher matrix of the 
power spectrum. For the spectrum plots in this section, the errors are based on the diagonal 
elements of the covariance matrix which is evaluated in Appendix iBl 

Figure [16] shows the effect that correlated noise has on the low i errors in the EE and BB 
spectra. The curves show the diagonal elements of the inverse Fisher matrix (the Ce errors) 
computed in two ways: (1) assuming the noise is uncorrelated in pixel space and described 
by N t> s (red) and (2) assuming it is correlated and correctly described by N -1 (black). The 
smooth rise in both curves toward low I is due to the effects of 1/f noise and is most 
pronounced in the W4 DA, which has the highest 1/f noise. The structure in the black 
trace is primarily due to the scan strategy. Note in particular, that we expect relatively 
larger error bars on £ = 2, 5, 7 in EE and on I = 3 in BB. We caution those analyzing 
maps that to obtain accurate results, the N _1 weighting must be used when working with 
the I < 23 power spectra. For the Monte Carlo Markov Chains (MCMC) and cosmological 
parameter evaluation, we do not use the power spectrum but find the exact likelihood of 
the temperature and polarization maps given the cosmological parameters (Appendix [D] & 



Hinshaw et all 120061 ). 



For both r4 and r9 maps there are 15 MASTER cross power spectra (see Table |5]). 
For the full three-year result, we form 5^f,- =1 ui x v j/6 omitting the ul x ul, u2 x u2, and 
u3 x u3 auto power spectra. In this expression, u and v denote the frequency band (K- 
W) and i and j denote the year. The noise per I in the limit of no celestial signal, Ni, is 
determined from analytical models that are informed by full simulations for r9 (including 
1/f noise), and from the full map solution for r4. 
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5.1. Power Spectrum of Foreground Emission Outside the P06 Mask. 

Figure [IT] shows the EE and BB power spectra for the region outside the P06 mask, 
74.3% of the sky, before any cleaning. The 15 cross spectra have been frequency averaged 
into four groups (Table [5]) by weighting with the diagonal elements of the covariance matrix. 
Data are similarly binned over the indicated ranges of I. It is clear that even on the cut sky 
the foreground emission is non negligible. In K band, we find £(£+1)C^ <2 _q > /2'k = 66 (fiK) 2 
and £(£+1)C^ <2 _ 6> /2tt = 48 (fiK) 2 , where £ =< 2 — 6 > denotes the weighted average over 
multipoles two through six. The emission drops by roughly a factor of 200 in Gt by 61 GHz 
resulting in < 0.3 (fiK) 2 for both EE and BB. There is a "window" between £ = 4 and 
£ = 8 in the EE where the emission is comparable to, though larger than, the detector noise. 
Unfortunately, BB foreground emission dominates a fiducial r = 0.3, r = 0.09 model by 
roughly an order of magnitude at £ < 30. In general, the power spectrum of the foreground 
emission scales approximately as £ _1//2 in £[£ + l)C?. 

Figure [TS] shows the power spectra as a function of frequency for a few £ bands. The 
spectrum of the emission follows that of synchrotron with T oc with (3 S = —2.9 for 
both EE and BE0. There is some evidence for another component at v > 60 as seen in 
the flattening of the EE £ = 2 term. We interpret this as due to dust emission. In the 
foreground model, we explicitly fit to a dust template and detect polarized dust emission. 
However, there is not yet a sufficiently high signal to noise ratio to strongly constrain the 
dust index or amplitude outside the P06 mask. 

A simple parameterization of the foreground emission outside the P06 mask region is 
given by 

£(£ + l)C/ or 727T = (B 8 (v/65) 20 ° + B d (is/Q5) 2f3d )£ m . (25) 

We have introduced the notation B xx = £(£ + 1)C XX /2tc to simplify the expression. The 
"d" and "s" subscripts stand for "dust" and "synchrotron." From an unweighted fit to all 
the raw £ < 100 data with the dust index fixed at (3d — 1-5, we find for EE B s = 0.36 (/iK) 2 , 
(3 S = -3.0, B d = 1.0 (/xK) 2 and m = -0.6; and for BB B s = 0.30 (/iK) 2 , (3 S = -2.8, 
Bd = 0.50 (yuK) 2 and m = —0.6. This model is given as an approximate guide. Its £ 
dependence is shown in Figure [T7J for v = 65 GHz and its frequency dependence is shown 
in Figure [T8l for BB £ = 2. One can see that this scaling model picks up the general trends 
but not the details of the foreground emission. For example, it ignores correlations between 
dust and synchrotron emission. It predicts an average foreground emission of ~ 1 (/iK) 2 at 



11 The fits to the power law index were done with the power spectra in CMB temperature units. The ratio 
of these and the indicies corresponding to antenna are approximately 1.03 and 0.99 at 90 GHz (where the 
difference is largest) for dust and synchrotron respectively. The difference is negligible. 
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30 GHz and I = 300. iLeitch et al.l (120051 ) give an upper bound of ~ 1 (/iK) 2 for synchrotron 
emission in this range. As DASI observes a relatively synchrotron-free region and at is 
beyond where this simple parametrization can b e tested, there is not a conflict with their 



results. The same is true for the CBI experiment (IReadhead et al.ll2004l ) which also observed 
at 30 GHz but at a predominantly higher I and in a predetermined clean region of sky. The 
dust amplitude in the model is especially uncertain. Depending on the region of sky within 
the P06 cut, and the £ of interest, it may be an order of magnitude off. 

For a more complete model of the power spectra of foreground emission, one must take 
into account the correlations or anticorrelations between various foreground components and 
between the foreground components and the CMB. For example, a reasonable fit to the 1 = 1 
EE spectrum, which is dominated by foreground emission, is given by 

B EE {v) = a.iyxvif' + p sd a s a d {^u^ + v{ A v^) + a d {v x v 2 f d (26) 

where p sd is the dust synchrotron correlation coefficient, V\ and are the frequencies of the 
two spectra that are correlated, the [3 d and (3 S are the dust and synchrotron spectral indices, 
and v = ^Jv\V2- This fit is shown in Figure [18] After normalizing the frequency to 65 GHz, 
the following coefficients were found to reasonably represent the data: a s = 0.64, (3 S = —2.9, 
a d = 0.65, f3 d = 1.5, and p sd = 0.46. In order to produce the KV, KW, and KaW features, 
there must be significant correlations between dust and synchrotron emission. For the I = 4 
EE spectrum a similar expression fits the data if p sd is negative. 

Some care is needed in interpreting the statistical significance of power spectra that 
include foreground emission and a cut sky. The lack of statistical isotropy of the foreground 
emission means that it must be treated separately from the CMB when assessing the net 
noise. In the presence of foregrounds, the random uncertainty becomes 

AC '-WTM- y lN ' +2N '^ ] (27) 

where Ti is the foreground emission at each £. We plot only the first term in Figures [17] 
& [18] to indicate the size of the statistical error. Additionally, with the sky cut there is a 
noise-foreground coupling between Nf' B and Ff^, and between N t ' and ■ This is 
analogous to the noise coupling shown in Appendix O 



5.2. Power Spectrum of Foreground-Cleaned Maps Outside the P06 Mask 



We next discuss the power spectrum after removing the foreground emission from the 
maps. Cleaning foregrounds not only changes the mean of Ci, but it reduces ACi because 
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of the couplings. The choice of model makes little difference to the conclusions. For all the 
following we have subtracted the best fit KD3Pol Q and U templates from the Ka through 
W maps (both r4 and r9 versions) as described in Section HI Table shows the EE 1 = 2 
and BB i = 5, the multipoles with the largest foreground contributions, for both before and 
after the subtraction. Where the foreground signal is dominant, the subtraction can reduce 
its level by a factor of 6-10 in temperature. 

When we fit and subtract the foreground templates, we use essentially all of the available 
data on polarized foreground emission. The error bar on the power spectrum of the cleaned 
maps is dilated in the cleaning process as discussed above. We do not include an additional 
error for systematic uncertainty in the model. Rather, by comparing spectra of pre-cleaned 
to cleaned maps, we estimate that the model removes at least 85% of the synchrotron. 
This is demonstrated, for example, in the KKa and KaKa combinations for i = 2 EE in 
Table [5j in the subtraction shown in Figure [T5J, and to a lesser degree by the null EB and 
BB power spectra. We also note that to a good approximation foreground emission adds 
only in quadrature to CMB emission. 

Figure dH] shows the power spectra of the foreground cleaned maps as a function of 
frequency for i = 2 — 9. It also shows what we estimate to be the maximum levels of residual 
foreground contamination in the power spectrum. In the figure, we plot the synchrotron 
spectrum scaled to 0.15 of the pre-cleaned Ka band value (in temperature). This shows that 
there is negligible residual synchrotron from 40 to 60 GHz with the possible exception of 
i = 2 at 40 GHz. Given the size of the 1 = 2 error bar, this potential contribution to the 
determination of the optical depth is negligible as discussed in Section §6.11 Constraining 
the residual dust contamination is more difficult. In Figure [191 we a l so show the MEM 
temperature dust model scaled by 5%, a typical dust polarization value. A similarly scaled 
FDS model is almost identical. This shows that even if we did not model and subtract dust, 
the contamination from it would not be large in Q and V bands. A more detailed model 
might have to take into account the possibility that the electrons and dust grains are in 
regions at different line of sight distances with different magnetic fields or that variations in 
the magnetic field could alias power from low multipoles to higher ones. 

The cross power spectra of the cleaned maps are combined by frequency band for testing 
cosmological models. The 10 cross spectra (since K-band is used in the model, there is no K- 
band cleaned spectrum) are assessed £ by t with a least squares fit to a flat line in Figure [191 
The results are shown for the QQ+QV+VV (denoted "QV combination") and QV+KaV 
combinations in EE, we find 0.1 < PTE < 1 for all i < 16, where PTE is for "Probability 
to Exceed" and is the probability that a random variable drawn from the same distribution 
exceeds the measured value of \ 2 - When W band is added to the mix, we find PTE < 0.03 
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for £ = 5,7, 9, though all other values of £ give reasonable values. For BB, all frequency 
combinations yield reasonable PTEs for all £. Thus, there is a residual signal in our power 
spectra that we do not yet understand. It is evident in W band in EE at I = 7 and to a 
lesser degree at £ = 5 and £ = 9. We see no clear evidence of it anywhere else. 

5.3. Null Tests and Systematic Checks Outside the P06 Mask 

Null tests are critical for assessing the quality of the data. We have examined the data 
in a wide variety of ways based on differencing assembly, frequency band, £ range, and year. 
We present selected, though typical, results in the following. A particularly important test 
is the null measurement of the BB, TB, and EB signals as shown in Tableland Figure [201 
These data combinations are derived from the same processing as the EE, TE, and TT 
combinations, where a signal is detected. Thus, the null result highlights the stability of the 
WMAP data, the mapmaking, the foreground cleaning, and the power spectrum estimation. 

The power spectrum of the difference of the individual yearly maps is another significant 
test. Table [7] shows the results for all the yearly differences for i = 2 — 16, the critical region of 
£-space for the cosmological analysis. We have also used the (ul x vl + u2 x v2 — ul x v2 — u2 x vl)/4 
cross spectra to similar effect. Here again the u and v denote different frequency bands. This 
combination is equivalent to forming the power spectrum of the difference between year one 
and year two maps. In principle it does not contain any signal. The cross-spectrum method 
treats the noise in a slightly different way from the straight map method, where one must 
use the error bars from one of the maps. It has been checked with simulations. Similar 
combinations are used for the other years. 

Using a variant of cross-spectrum method, we have also tested combinations of DAs for 
multiple ranges in £ within each frequency band. For all null tests, we find the expected 
null measurements, apart from the previously mentioned residuals at £ = 5&7 in W band. 
Table Ogives the reduced x 2 for all combinations of T, E, and B data for a number of data 
combinations. 

From Figure [191 h is clear that the large signal in W band is not residual dust contam- 
ination because the dust would not fit measurements in VW. Additionally, if one assumes 
that the polarized emission at a particular £ is a fraction times the intensity at the same £, it 
would require > 40% dust polarization, which is unreasonable. Though this simple picture 
does not take into account the aliasing of intensity from a lower £, we do not observe a similar 
effect with the synchrotron emission, which in the simplest case is polarized by the same 
magnetic fields. The W-band EE £ = 7 value is essentially unchanged by cleaning, removing 
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a 10° radius around the Galactic caps, or by additionally masking ±10° in the ecliptic plane. 

A number of tests have been done to identify this artifact of the data. We are not yet 
certain if it is due to an ersatz signal or an incorrect noise term. The error bars on the 
individual year differences are too large to clearly see if the effect is the same from year 
to year. Simulations show that 1/f alone cannot explain the signal. The scan pattern in 
combination with the change in polarization is directly related to the large error bars at 
I = 5, 7, 9 and is well understood. We have not identified a mechanism that leads to a 
further increase in these uncertainties. We know that different treatments of the noise, for 
example using N f, s weighting, decreases the magnitude of the discrepancy, though we are 
confident that the N -1 treatment of the pixel noise is the correct approach. The discrepancy 
can be made smaller by eliminating the Wl data simply because the error bars increase. The 
Wl radiometer has the lowe st noise but also the largest number of " glitches" (13, 4, 1 in 



years 1, 2, & 3 respectively, iLimon et al.l 120031 ; iHinshaw et al.ll2003bl ). However, since we 



could not identify any correlation between the glitch rate (assuming that unmasked glitches 
are responsible) and the magnitude of the signal, we do not have a basis for eliminating this 
channel. 

We believe there is an as yet unknown coupling in the W-band data that is driving the 
signal but more simulations and more sensitivity are needed to understand it. We cannot 
rule out similar lower-level problems in other bands, but we see no evidence of systematic 
effects in BB, EB, or other values of I and other frequencies in EE. To avoid biasing the result 
by this residual artifact which also possibly masks some unmodeled dust and synchrotron 
contamination, we limit the cosmological analysis to the QV combination. We also show 
that including W band EE does not alter our conclusions. 



5.4. Analysis of Foreground-Cleaned Power Spectra Outside the P06 Mask. 

A comparison of the raw spectra and foreground cleaned spectra is shown in Figure [2TJ 
We start with the weighted sum of the 8 cross spectra with v > 40 GHz (without KW). 
This is the upper-level line (green) in the figure. The individual maps are then cleaned and 
the power spectra remade and coadded. This is shown in violet. Similar comparisons are 
repeated for the QVW and QV combinations. A simple visual inspection shows that even 
at the is with the highest foreground contamination, the cleaning is effective. 

From the bottom left panel in Figure [21] one sees that there is a clear signal above 
the noise in EE at i < 7. For the QV combination, B%? <2 _ 6> = 0.086 ± 0.029 {fiK) 2 . 
The signal has persisted through a number of different analyses. We cannot rule out that 



-33- 



this signal might find explanation in an unmodeled foreground component; however, we 
find this explanation unlikely since the emission would have to be strikingly different from 
the measured spatial and frequency characteristics of the polarized foreground emission. 
Additionally, when different bands are coadded, the signal level is consistent: for QVW 
0.098 ± 0.022 (/jK) 2 and for all channels with v > 40 GHz except KW, 
0.095 ±0.019 (nK) 2 . We have searched for systematic effects in the EE £ = 2 — 8 
range and have not been able to identify any, other than the one discussed above. We can- 
not find a more plausible explanation than that the signal is in the sky. We are thus led to 
interpret it cosmologically. This is done in the next section. 

We show the EE signal for £ > 20 in Table M and in Figure [221 along with a comparison 



Y2EE 

u £=<2-6> 
K>EE 

u l=<2-6> 



to other recent mea surements (ILeitch et al.l 120051 ; ISievers et al.l 120051 ; iBarkats et al.l 12005 



Montroy et al.ll2005l ). Based on the best fit to the TT spectra, we produce a template for 



the predicted EE spectrum Cf E ' T and form: 



800 



(28) 



1=50 



where SC t = Cf E - A EE C EE ' T , A EE is the fit amplitude, and Qf e E is the diagonal Fisher 
matrix in Appendix I C. 3 1 Off diagonal elements in Qf E have a negligible effect on the results. 



The results of the fit are plotted in Figure [23] for various frequency combinations. We 
plot Ax 2 from the minimum value and find that A EE = 0.95 ± 0.35 for the pre-cleaned 
QVW combination, where the uncertainty is determined from the bounds at A% 2 = 1. The 
reduced \ 2 at the minima are 1.34, 1.34, 1.24, and 1.30 for QV, VW, QVW, and KaQVW 
combinations respectively, most likely indicating residual foreground contamination. At a 
relative amplitude of zero, A% 2 = 1.0, 3.3, 6.2, & 16 respectively for the same frequency 
combinations. It is clear that the noise is not yet low enough to use just QV as was done 
at low multipoles. In addition, cleaning the maps with the KD3Pol is problematic because 
the K-band window function is reduced to 0.1 by £ = 250. When the same code is used to 
analyze EB and BB data, the fitted amplitude is always consistent with zero. To summarize, 
the WMAP EE data are consistent with a model of adiabatic fluctuations based on the 
temperature maps at greater than the 2a level for the QVW and KaQVW combinations. 



Figure [24] (left panel) shows the TE spectrum for I < 16. We use V band for tempera- 
ture and the QV combination for polarization. Several aspects of the new processing led to 
increased errors an d a reduced \ow-£ signal estimate relative to the first-year result (Figure 8, 
Kogut et al.ll2003l ). These include: improvements in mapmaking and power spectrum esti- 
mation (especially accounting for correlated noise and applying N -1 weighting); limiting the 
bands to just Q and V instead of Ka-W; increasing the cut from KP0 to P06; and improve- 
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ments in foreground modeling, including a new estimate of dust polarization. Recall also 
that the first-year result was based on the combinations of Ka, Q, V, and W bands and did 
not include a dust polarization template in contrast to the new prescription. Furthermore, if 
the year-two data are processe d in th e same way as the first-year data, we obtain a spectrum 
similar to that in lKogut et al.l (120031 ) indicating that the major difference between first-year 
and three- year results rests on new knowledge of how to make and clean polarization maps. 
The new spectrum is fully consistent with the first-year results and prefers a model based 
just on TT and EE data to a null signal at the 2a level. However, the new spectrum is 
also consistent with the absence of a TE signal. Thus, it will take greater signal-to-noise to 
clearly identify the TE signal with our new analysis methods. 



Figure [2H (right panel) shows the TE signal ove r the full range in £. Other detections 
of TE at £ > 100 have been reported by D ASI (2.9a) jLeitch et a DlaOOpj), Boomerang (3.5a) 



Jpiacentini et alJbood ). and CBI (3.3a) Jsievers et al.l l2005h . The WMAP data have had 



foreground models subtracted from both the temperature and polarization maps prior to 
forming the cross correlation. The expected anticorrelation between the polarization and 
temperature is clearly evident. To quantify the consistency with the TT data we make a TE 
template based on the model fit to TT. Next, a fit is made to the TE data for 20 < £ < 500 
with the following: 

X 2 (A TE , M) = J2 SC e Qj £ ?6C £ , (29) 

w 



C, ' J (Af), Cf ' (A£) is the predicted power spectrum shifted by 



where 5C e = Cj E - A TE 
A£, A TE is the fit amplitude, and Qj E is the diagonal Fisher matrix in Appendix 10.31 Off 
di agonal elements i n Qj ( E have a negligible effect on the results. Similar 2D fits were done 



in 



Readhead et al.l (120041 ) . We show the combination that uses V and W bands for T and Q 
and V bands for E. The result, shown in Figure[23l is A TE = 0.93±0.12 and A£ = 0±8 with 



X 2 /V = 468/482 (PTE=0.66). Similar results are obtained with other band combinations. 
Thus the TE data are consistent with the TT data to within the limits of measurement. 



Figure [25] shows a summary of the various components of the CM B anisotropy. 



6. Cosmological Analysis 

The £ < 100 region of the CMB polarization spectra is rich with new tests of cosmology. 
The EE spectrum gives us a new measure of the optical depth. The same free electrons from 
reionization that lead to the £ < 10 EE signal act as test particles that scatter the quadrupo- 
lar temperature anisotropy produced by gravitational waves (tensor modes) originating at 
the birth of the universe. The scatter results in polarization B modes. Tensor modes also 
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affect the TT spectrum in this region. A combination of these and related observations leads 
to direct tests of models of inflation. 



The detection of the TE anticorrelation near 



30 is a fundamental measurement 



of the physics of the formation of cosmological perturbations (IPeiris et al.l 120031 ) . It re- 
quires some mechanism l ike inflation to produce and shows that superhorizon fluctuations 



must exist. iTurokl ( 119961 ) showed that with enough free parameters one could in principle 
make a model based o n pos t-inflation causal physics that reproduced the TT spectrum. 
Spergel fc Zaldarriagal (119971 ) show that the TE anticorrelation is characteristic of models 
with superhorizon fluctuations. The reason is that the anticorrelation is observed on angular 
scales larger than the acoustic horizon at decoupling. Thus, the observed velocity-density 
correlations implied by the TE data must have existed on scales larger than the horizon and 
were not produced by post-inflation causal processes. 

Although multiple distinct physical mechanisms affect the I < 100 spectra, their effect s 
can be disentangled through an analysis of the full data complement (ISpergel et al.l I2006T ) . 
The separation, though, is not perfect and there remain degeneracies. In particular, to some 
degree, the values of the scalar spectral index, n s , optical depth, and the tensor to scalar 
ratio, r, may be traded against each other, although far less than in the first-year WMAP 
results. As the data improve, or as more data sets are added, the degener acy is broken 



furth er. In the following we take a step back from the full MCMC analysis (ISpergel et al. 
20061 ) and estimate r and r from analyses of just the i < 10 polarization spectra. This 



approach aids our intuition in understanding what it is in the data that constrains the 
cosmological parameters. 



6.1. The Optical Depth of Reionization 

Our knowledge of the optical depth ripples through the assessment of all the cosmic 
parameters. Free electrons scatter the CMB photons thereby reducing the amplitude of the 
CMB spectrum. This in turn directly impacts the determination of other parameters. 

The distinctive signature of reionization is at i < 10 in EE. The only known contami- 
nation is from foreground emission which has been modeled and subtracted. The amplitude 
of the reionization signal is proportional to r in TE and is proportional to r 2 i n EE and BB. 



In the first year analysis, we imposed a prior that r < 0.3 (ISpergel et al.ll2003l ). Such a high 



value would produce a signal > 6 times the model in Figure [21] and is clearly inconsistent 
with the EE data. Thus this new analysis is a significant improvement over the previously 
assumed prior. 
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We assess r using three methods: (1) with template fits to the EE power spectra; (2) with 
an exact likelihood technique based directly on the maps as de scribed in AppendixlD| and (3) 



with a multiparameter MCMC fit to all the data as re ported in ISper gel et al.l (120061 ) . The first 



method is based directly on the MASTER spectrum (IHivon et al.ll2002l . and Appendix [B]) of 
EE data and serves as a simple check of the other two. Additionally, the simplicity allows 
us to examine the robustness of the EE and TE detections to cuts of the data. The second 
method is robust and takes into account the phases of the EE and TE signals. It is run either 
as a stand alone meth od, as reported here, or as part of the full MCMC chain as reported 
in ISpergel et al.l (120061 ) . The best estimate of the optical depth comes from the full chains. 



For the template fits, ACDM power spectra were generated for < r < 0.3, with the 
remaining parameters fixed to n s = 0.96, a;& = 0.0226, u m = 0.133, and h = 0.72. For each 
spectrum, the scalar amplitude A is fixed by requiring that Bj^ 200 = 5589 ^K 2 . We then 
form: 

C(f) = ; exp[- - x t I l )D- 1 (x l - ff)/2] (30) 

(27r)« v /det( J D) V 

where xi is the data as shown in Figure I2T1 xf 1 = Bf E (f) is the model ACDM spectrum, 

D ' = wniiw (B?E(f) + NF? (31) 

as in C.14, and N EE is the uncertainty shown in Figure [21] and is derived from the MASTER 
spectrum determination. We use the symbol f in this context because the likelihood function 
we obtain is not the full likelihood for r. Uncertainties in other parameters, especially n s , 
have been ignored and the Ce distribution is taken to be Gaussian. Thus C(f) does not give 
a good estimate of the uncertainty. Its primary use is as a simple parametrization of the 
data. We call this method "simple tau." Table [9] shows that simple tau is stable with data 
selection. One can also see that if the QQ component is removed from the QV combination, 
t increases slightly. This is another indication that foreground emission is not biasing the 
result. Additionally, one can see that removing i = 5, 7 for all band combinations does not 
greatly affect r. 

The optimal method for computing the optical depth is with the exact likelihood (as in 
Appendix [D]). The primary benefits are: it makes no assumptions about the distribution of 
Cg at each I but does assume that the polarization signal and noise in the maps are normally 
distributed; it works directly in pixel space, taking advantage of the phase relations between 
the T and E modes both together and separately; and it is unbiased. The only disadvantages 
are that it is computationally intensive and that it is not easy to excise individual values of 
t such as t = 5, 7. 

In the exact likelihood method we take into account errors in our foreground model by 
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marginalizing the likelihood function given in (Appendix [D]) over the errors in the fitting 
coefficient for synchrotron emission, a s . We ignore foreground errors in dust emission, as 
polarized dust emission is negligible in any of combinations of frequencies reported in Table 9. 
When errors in a s are Gaussian, the marginalization simply yields an additional term in the 
covariance matrix, 

Cij = + > Cij = + + Fij, (32) 

where 

Fij = alfifj, (33) 

and a 2 s = {a 2 s ) — (a s ) 2 and /, is a template map of polarized foreground (i.e., Q and U 
maps in K band). Here, the mean values, (a s ), are given in the second column of Table 3. 
We find a s = 0.007 in Q and V band. To be conservative, we adopt a s = 0.01 as our 
foreground error, which is the 2-a bound on the foreground error in QV combination. As the 
foreground marginalization yields a new positive term in the covariance matrix, a fraction 
of the signal that was attributed to CMB before is now attributed to foreground, when the 
spatial distribution of the signal is the same as that of K-band maps. The values of r with 
the foreground marginalization are tabulated Table 9. The marginalization reduces r by 
0.0017 in QV. The largest effect is seen in Q band, for which r drops by 0.0027. Thus, the 
foreground error does not significantly affect our determination of r. 

Table [9] shows that similar vales of r are obtained for a wide variety of band combina- 
tions. This is another indication that foreground emission is not significant. We conserva- 
tively select the QV combination. Table M also compares the exact likelihood for the EE QV 
combination to the simple tau method. One can see that simple tau is slightly biased high 
when compared to the exact likelihood and underestimates the likelihood at r = 0. One 
source of the bias is the assumption of a Gaussian likelihood. Nevertheless, it is reassuring 
that a variety of combinations of data give consistent values of r. 

The values given here are just for the EE and TE data considered alone, with the 
first peak TT ampli tude fixed. When the exact likelihood is used in the full MCMC analysis 



( ISpergel et al.ll2006l ) yielding the best estimate, we find r = 0.089 ±0.030, slightly lower than 
the values reported here but with the same uncertainty, indicating that the simple analysis 
has exhausted most of the information on the optical depth contained in the polarization 
data. 



As discussed in (ISpergel et al.l 120061 ) there is a degeneracy between the scalar spectral 
index, n s and r. If we had instead selected the K-band directions for the dust polarization 
template, we would have found r = 0.107 and an increase in n s of 0.004. A similar shift 
would have been found using the KaQVW combination shown in Table M and Figure [2HJ 
This is another indication of the relative insensitivity of the results to the cleaning method. 
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6.2. Gravitational Waves 



The Cf B spectrum directly probes the primordial gravitational wave background pro- 



was proposed by 


Starobinskv 


dl979 


). 


Liddle & Lvth 


( 


200ol): 


Dodelsonl ( 


2003) 



ations both contribute to the TT and EE spectra, only tensors produce B modes. Inflation 
models generally predict similar scalar spectra, bu t differ in their prediction of the tensor 



comp onent. For example, ekpyrotic/cyclic models (IKhoury et al.ll2002l ; ISteinhardt fc Turok 



20021 ) predict no observable tensor modes. 



The tensor contribution is quantified with t he tensor to scalar ratio r. We follow the 
convention in the CAMB co de Lewis et al.lboool. Vers i on, June 200 4), in CMBFAST v4.5.1 
JSeliak fc Zaldarriagalfl9~9~6r i and in lPeiris et al.l J2OO3I ): IVerde et al.l (bo03h : 



Aj(fco) 
A^ ) 



(34) 



Here, A^ and A^ are the variance due to scalar and tensor modes respectively. They are 
defined through 



(II 2 ) 

and (h% im h prim < ij ) 



dk 
~k 



(35) 
(36) 



where h%p m is the primordial tensor metric perturbation in real spac e that was generated 
during inflation and stretched to outside the horizorl 1 " 2 !. Peiris et al. ( 2003 ) shows the k- 
dependence of these expressions. 

The expression for r is evaluated at k = 0.002 Mpc -1 correspond ing to I ~ rink = 3 
with the distance to the decoupling surface rj « 14, 400 Mpc. Following IVerde et al.l (120031 ). 
we use A^(fco) = 2 x 10 4 7r 2 A(/c )/9T 2 w 2.95 x lQ- g A(k ) with T in microk elvins. Some of 
the s i mple models of in flation in a ACDM cosmology predict r ~ 0.3 (e.g., iLiddle fc Lyth 



20001 ; iBoyle et al.l 120051 ) . For example, near this range inflationary model s with a ma ssive 
scalar field, V{(p) = m 2 <p 2 /2, predict r = 8/N e = 4(1 - n s ) = 0.13 - 0.16 JLinde!ll983l ) and 
models with a self coupling, V(<f>) = A0 4 /4, predict r = 16/A" e = 16(1 — n s )/3 = 0.27 — 0.32 
for N e = 60 — 50. Here, N e is the number of e-foldings before the end of inflation. However, 



12 Note that our convention yields r = 16e for slow-roll inflationary models with a single scalar field. Here, 
e is the slow roll parameter related to the square of the slope of the inflaton potential. 
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some variants produce r > 0.32 (e.g. iMukhanov fc Vikmanl 120051 ) while many other have 
r < 0.1. 



For the best fit WMAP-only ACDM plus tensor model, the optical depth is r = 0.091. 

If we add to this model a tensor component with r = 0.3, then Bf B <2 _ &> = 0.001 /iK 2 . 

3 _ 
<2-6> — 



A simple average of the Cf B data gives Bf B <2 _ 6> = -0.044 ± 0.030 /iK 2 , B BB 



-0.018 ± 0.023 /iK 2 , Bf B <2 _ 6> = 0.003 ± 0.020 /iK 2 , for QV, QVW, and v > 40 GHz (no 
KW) combinations respectively. To detect a signal at the upper range of the predictions 
would require maps with « 5 times smaller error bars. 

We constrain r by directly fitting a template of Cf B to the BB data. With the above 
definition, r directly scales the C BB power spectrum. Additionally, the amp litude of C B B 



for I < 16 scales as r 2 . We set the template to be the standard ACDM model (jSpergel et al. 



20061 ) and use the single field inflation consistency relation, n t = — r/8, to fix the tensor 
spectral index. We assume the spectral index does not run and set n s = 0.96. We distinguish 
the r in the template fit by the f notation. The sum is over 2 < £ < 11. 

The results of the fit are plotted in Figure [23 When we consider just the limit on f 
from the polarization spectra, ignoring the tensor contribution to TT, we find f < 2.2 (95% 
CL) after marginalizing over f. It is clear that the BB spectrum is not driving the limit on 
r. After including the TT data, the limit drops to f < 0.27 (95% CL). This shows that the 
TT data in combination with the limits on r from EE and TE are leading to the limit on r. 
The full MCMC analysis gives r < 0.55 (95% CL) with just the WMAP data. The increase 
in the error over the simple method given above is the result of the marginalization over the 
other parameters, particularly n s . Additionally, when n s is allowed to depend on k, the error 
in r increases dramatically, allowing r < 1.3 (95% CL). 

We can relate r to the curren t energy density in primordial gravitational radiation 
rauss fc Whitelll992l ; |Peirisll2003h . 

^ow = ^iJ yA 2 (fc)T 2 (M), (37) 

where 7] is conformal time and the transfer function, T(k,r]), is given in Equation lE18l 
The approximation given in Equation IE31I evaluated for A = 0.838 and f < 2.2 yields 
tt GW < 9.6 x 10" 12 (95% CL) and for r < 0.55, Vt GW < 2.0 x 10~ 12 (95% CL). 



7. Discussion and Conclusions 



WMAP detects significant levels of polarized foreground emission over much of the 
sky. The minimum in contamination is near 60 GHz outside the P06 mask. To detect the 
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polarization in the CMB at I < 10 a model of the foreground emission must be subtracted 
from the data. This situation differs from that of the analysis of the temperature anisotropy 
for which the foreground emission may be simply masked as a first approximation^ 

WMAP has detected the primary temperature anisotropy, the temperature polarization 
cross correlation, and the E-mode polarization of the CMB. We detect the optical depth with 
r = 0.089 ± 0.030 in a full fit to all WMAP data. This result is supported by stand-alone 
analyses of the polarization data. Using primarily the TT spectrum, along with the optical 
depth established with the TE and EE spectra, the tensor to scalar ratio is limited to 
ro.002 < 0.55 (95% CL). When the large scale structure power spectrum is added to the 



mix (jSpergel et al.ll2006l ). the limit tightens to ro.002 < 0.28 (95% CL). These values are 



approaching the predictions of the simplest inflation models. 

A clear detection of the B modes at £ < 100 would give a direct handle on the physics 
of the early universe at energy scales of 10 15 — 10 16 GeV. This paper shows that care will 
be required to unambiguously separate the intrinsic signal from the foreground emission. 
However, the BB spectrum is particularly clean in WMAP and, at least for £ = 2,3, the 
foreground contamination is relatively low. In the noise dominated regime, the error bar on 
Cf B decreases in proportion to time. Continued WMAP operations combined with other 
experimental efforts are nearing a range of great interest. 

These new results involve a complete reevaluation of all the components of our previous 
analyses, from the beams and gain models through to the mapmaking and foreground mod- 
eling. The data and most of the derived data products are available through the LAMBDA 



website, http://lambda.gsfc.nasa.gov/ WMAP continues to operate nominally. In the fu- 
ture we will address a number of the open issues raised above. In particular, we can anticipate 
a better understanding of systematic errors and foreground emission, and therefore improved 
constraints on r and r. It is remarkable that our understanding of the cosmos has reached 
the point where we have begun to quantitatively distinguish between different models of the 
birth of the universe. 



8. Acknowledgments 

We wish to thank Uros Seljak, Suzanne Staggs, and Paul Steinhardt for enlightening 
conversations on polarization and inflation, and Jim Peebles for his explanation of reioniza- 
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((Bennett et al.ll2003bl ). 



-41 - 



tion. We thank Bruce Draine for helpful discussions on dust polarization and Rainer Beck 
for many suggestions on modeling the Galactic magnetic field. We thank Andre Waelkens 
and Torsten Enfilin for pointing out a mistake in the original calculation of g sync . (The fix 
affected Figures 8 and 9 but not the conclusions.) We thank Bruce Winstein for a detailed 
set of comments on the submitted draft. EK acknowledges support from an Alfred P. Sloan 
Research Fellowship. HVP wishes to acknowledge useful discussions with R. Easther, S. 
Larson, D. Mortlock, and A. Lewis. The WMAP mission is made possible by the support 
of the Office of Space Sciences at NASA Headquarters and by the hard and capable work 
of scores of scientists, engineers, technicians, machinists, data analysts, budget analysts, 
managers, administrative staff, and reviewers. HVP is supported by NASA through Hubble 
Fellowship grant #HF-01177.01-A awarded by the Space Telescope Science Institute which 
is operated by the Association of Universities for Research in Astronomy, Inc., for NASA 
under contract NAS 5-26555. This research has made use of NASA's Astrophysics Data 
System Bibliographic Services, the HEALPix software, CAMB software, and the CMBFAST 
software. This research was additionally supported by NASA LTSA03-000-0090, NASA 
ATPNNG04GK55G, and NASA ADP03-0000-092 awards. 



A. Radiometer Model 

In this section we develop a simple model for the WMAP instrument using Jones matri- 



ces (Jones 



1996 



1941 



Mo ntgomery et al.lll948l ; lBlumlll959l ; lFarislll967l ; ISault et al.lll996l ; iTinbergen 



Hu et al.l 120031 ). In the following we assume that all circuit elements are matched and 



ignore additive noise terms. 

The Jones matrix J models the instrumental response to polarization, 



E 



out 



JEu 



(Al) 



linearly relating the output electric field to the input. WMAP is a differential instrument, 
so the input radiation vector E [n has four elements, (Eg-, Ey, E^ , Ey), corresponding to the 
electric field seen by the A- and B-side feed pair. The outputs E ont are the inputs to the 
detectors. 



The first link in the chain is to model the optics, feeds, and orthomode transducers 
(OMTs). We consider them as a single unit, because ascribing effects to the individual 
components is difficult and not well defined in terms of observations. We include two effects, 
loss imbalance and polarization leakage: 



rA,B 
'OFO 



-rA,BjA,B 

loss crosspol 



(A2) 
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y 



J loss - ( q L A,B ) ( A3 ) 



Xt' B e^ B 



J crosspol - ( _X^' B e' iY ^ B 1 1 ^ 

Here L£'y is the loss for the particular polarization and X^ quantifies the level of cross- 
polarization (or polarization isolation) leakage, which we model as a small rotation error. 
The matrix J^^po! ^ s the ^ TS ^ term in the expansion of a general unitary matrix but is not 
unitary itself. The subscripts "1" and "2" refer to the two orthogonally polarized radiometers 
which are differenced to form Ap. The matrix J^fspoi * s the ^ TS ^ term in the expansion 
of a general unitary matrix, though it is not unitary itself. The cross-polarization terms 
are allowed to have arbitrary phases Y 1 2 ' . It is possible for cross polarization to produce 
circular polarization but WMAPcannot detect it in Aj or Ap. While in general there are 
four loss terms, two of them are calibrated out. The two that remain are the radiometer loss 



imbalances, x imj i and x- im ^- IJarosik et al.l (j2003bl . Table 3) measured the loss imbalances by 



fitting the response to the common mode CMB dipole signal, and found them to be < 1%. 
The mean imbalance, x im = (x imj i + x imj2 )/2, is corrected for by the map-ma king algorithm, 



while the "imbalance in the imbalance", Sx im = (x imj i — x im]2 )/2, is not (IHinshaw et al. 
2003al . §C3). To connect the different notations, = Li(l + 2i m ,i), Ly = L 2 (l + Xi mj 2 



Z/f = L 2 (l — x im 2 ), and Ly = Li(l — x imi i). The L x and L 2 are calibrated out. 



The next step is to model the radiometers. They are described in detail in lJarosik et al. 



(j2003al ). so we simply present the Jones representation of the radiometer and refer the reader 



to the paper for more details. 

Jradiometer JwarmT Jswitch Jamp JcoldT ( A ^) 

1 



' coldT 



1 switch 



' warmT 
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V2 ( 
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(A6) 



1 

J amP = ( n ) (^7) 

(A8) 
(A9) 



1 

-1 



Here, g s , gd are the amplifier gains in the two legs of the radiometer, and is the instantaneous 
phase of the (unjammed) phase switch. We have lumped the warm and cold amplifiers 
together. 
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^ bandpass Jradiometer-^- connect JqFO 



(AlO) 
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connect 



' radiometer 



' bandpass 




(All) 



(A12) 



(A13) 



The detector outputs in counts (013,014,023,024) are the diagonal elements of P out = 
(E out El ut ), multiplied by the responsivities (S13, S14, S23, S24). Jofo is a 4x4 matrix with 
Jq FO (Equation IA2j) filling the upper left 2x2 entries and Jqfo filling the lower right 2x2 
entries. 

Pout = JDAPinJk (A14) 

P A 

in 





(A15) 



T x + Q x 
U x + iV x 



U x ■ 

rpX 



iV x 

Q x 



(Ai6) 



In this expression, Stokes Q, U, and V refer to the quantities measured in the radiometer 
reference frame; we drop the "Rad" notation used in $3] for notational convenience. Before 
the outputs are recorded they are demodulated in phase with the phase switch. We model 
this process as 

h + vr - 5t)] (A17) 



1 

2" 



Cij{ 



where (pi is the phase difference between the two radiometer legs, and 5i is the error between 
the two switch states. 



Since the input radiation is incoherent, 



out 



duo 



OP 



out 



duo 



{Mi 
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Since Jbandpass is the only frequency dependent component in the model, we make the sub- 
stitution flj(uj) — > ffj, where 

H = J du${u). (A19) 

The calibrated detector outputs are dij = Cijj G^, where is a gain for the temperature 
difference, 

G^ = ^LigisgidftjSij cos(5j/2) cos(0j - 8i/2) (1 - e#). (A20) 

Here e»j is the calibration uncertainty. 

The radiometer signal channels are AT = (d$ 3 — da)/ 2, from which are formed the 
temperature and polarization signal channels AT/, ATp. Then to first order in the systematic 
uncertainties, 

AT/ = 2 5x im Q- + e~ Q + + (l + e + ) T~ + 2 x im T+ + Z A U A - Z B U B (A21) 
ATp = 2x im Q~ + (l + e + ) Q + + e~T- + 2 5x im T + + Z A U A + Z B U B (A22) 

Here T± = T A ± T B , {Q ± ,U ± ,L ± ,} are similarly defined, Z±' B = X A ' B cos(Y A ' B ) ± 
cos(F 2 A '' B ) encodes the influence of the crosspol effects, and e ± = ((ei 3 + e 14 ) ± (e 2 3 + 
£24))/4. The dominant ATp component is Q + , not Q~ , because Q A — > —Q B ,Q B —> —Q A 
when the spacecraft rotates 180°. In the limit of no loss imbalance or calibration error, and 
similar cross polarization for all components, ATp = Q + + 2Xcos(Y~)£/ + . 



B. Estimation of the Polarization Power Spectra 



The WMAP polarizat ion power spectra at / < 32 incorporate an extension of the MAS- 
TER quadratic estimator (IHivon et al.l l20Q2f ) . which is used to account for mode coupling. 
The original method assumes that observations of every point on the sky give statistically 
independent noise. However, WMAP has a significant component of the noise that is cor- 
related between pointings due to its scan pattern and the 1/ / noise, and thus the method 
needs to be modified as described here to accommodate a full covariance matrix. The most 
conspicuous mathematical feature of the original method is Wigner 3-j symbols, whereas in 
the extended method, these objects are not used. For more details of the origina l method, 
as well as the application to polarization, see Appendix A of iKogut et al.l (120031 ). together 
with the references therein. 
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B.l. Extended MASTER Algorithm for Temperature Power Spectrum 

The original method is derived by modeling the sky brightness as a continuous function 
of pointing. For example, the observed cut-sky spherical harmonic coefficients for Stokes /, 
denoted as T[ m , are defined as follows: 



Here, n is the unit vector of the pointing, w(n) is the weighting function, T(n) is the sky 
brightness, and Yi m (n) is a spherical harmonic basis function. Expanding T(n) and w(n) in 
spherical harmonics gives a series. Each term of the series includes an integral of a product 
of three spherical harmonic basis functions: 



These distinctive integrals are what give rise to the 3-j symbols. The orthogonality relations 
of 3-j symbols eliminate many terms in the expression for the observed power spectrum. 

When there is noise covariance, the weight is a function of two pointings rather than 
just one, and the 3-j symbols are not used. This case is most easily treated by modeling 
the sky as a set of discrete pixels. The goal of the derivation is to form a mode-coupling 
matrix M^, Y ' X Y , where XY and X'Y' are each chosen from the nine correlations TT, TE, 
TB, ET, EE, EB, BT, BE, and BB. In order to introduce the formalism, we first discuss 
the TT correlation, which is the simplest. Because there is no coupling between TT and 

TT TT 

the other eight correlations, only M n , ' needs to be considered. We note here that we do 
not actually use this formalism for TT but only for the others, as the temperature power 
spectrum at low-/ is dominated by the signal and the noise correlation is not important. We 
use TT here to illustrate the main point of the method. The extension to the polarization 
power spectra that follows TT (§ IB.2j) is what we use for the actual analysis. 

The weighting is computed initially as the inverse of the covariance matrix of the pixels. 
The sky cut is expressed by setting the appropriate rows and columns to a very large number 
in the noise covariance matrix before inverting it [Eq. (ID7I) ]. We call the resulting weight 
matrix W. Further, let Yi m ^ be a matrix containing (appropriately normalized) values of a 
spherical harmonic basis function evaluated at each pixel, p. index Im. The number of rows 
of Yim,p is n p , which is the number of pixels in each sky map. The observed Stokes / sky 
map is T p . 

In this notation, the observed spherical harmonic coefficients are expressed as 






dn Yi 
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If the matrix W is diagonal, this expression is simply the discrete version of Eq. IB1I above. 
Expanding T p i in spherical harmonics gives 



Tim — 



I'm' 



I'm' ,p' 



VP 



Tllrnf- 



This expression suggests the utility of defining 

Zl m ,l'm' = ^ ^ Ylm,pWpp'Yl'm' ,f 



(B2) 



so that 



VP 



Tim — ^ Zi rri ,l'm'Tl 
I'm' 



I'm' 



The value of the observed power spectrum at I is expressed as follows: 

(2J + l)C, = J2 f lmTlm 
m 

= 22 ^2 5Z^' m ' / " m "^" m "''*^ /m ' / ' m '"^ 



I'm' 



(B3) 



m Z"m" Z'm' 



In order to get the true, underlying CMB power spectrum into the equation, the next step 
is to take the expectation of Eq. 



(B4) 



(21 + 1) (<5j) - /2 z2 Z lm,l"m" Z lm,l'm' (^l*m"^'m') 

l"m" I'm' m 

= ^2 z2 z2 Z lm,l"m" Z lm,l'm' (Cl')6l'l"8 m >, 
l"m" I'm' m 

= ^2 I /j Z lmMm' Z lm,l'm' 1 (Q') ■ 
I' \mm' / 

Therefore, we obtain the unbiased estimator of the underlying power spectrum as 

^ = E( M ~%^> ( B5 ) 

where 



Mm 



21 + 



— ^2 \ Z lm,l'm' 



(B6) 



In order to apply this method to cross-correlations between DAs, one of the Z matrices in 
Eq. IB6I is computed from the noise matrix of the first DA, and the other from that of the 
second DA. 
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B.2. Extended MASTER Algorithm for Polarization Power Spectra 

The same formalism accommodates polarization. In what follows, uppercase X or Y 
indicates one of the three harmonic transforms T, E, or B, and lowercase a or b denotes 
the Stokes parameter label /, Q, or U . The following substitutions are made in the above 
derivation: 

Wppi — > W( aj) )( /p/) (B7) 

Ylm,p T(x«m)(ap) (B8) 

where the non-zero elements of T are 

Y(Tlm)(Ip) — Ylm,p (B9) 
T(_Bim)(Qp) = — g (+2^m,P + -2^m,p) (BIO) 

T(B/ m )(Q p ) = — - (+2^m,p — -2^m,p) (Bll) 
Y(Elm)(Up) — {Blm)(Qp) (B12) 

"f(Blm)(Up) — Y(Elm)(Qp) (B13) 



±2^m,p are spin-2 spherical harmonics in the same matrix form as Y\m, 



v 



For each pair of DAs, a Z matrix is computed by analogy with Eq. IB2I The derivation 
follows the general steps above. The analog of Eq. IB5I is 



Cr= £ (M- 1 )^' ~cf Y \ (B14) 



A/f XY,X'Y' _ 1 yXX'* rvYY 1 

1V1 W ~ 21 _|_ I / j ^lm,l'm lZj lm,l'm'i 



,X\ ... ^ ' ^r-l\ xvx ' ] ' r.A'V 
X'Y'l 

where 

M* Y > X ' Y ' = 

21 + 

where 



Z lL X l'm' = 22 ^\xi m ){ap)W(ap){a'p>)^(Xn> m >)(a>p>) ( B 15) 
ap,a'p' 

For each DA pair, the 81 coupling submatrices M XY ' X Y are combined in a grand 
coupling matrix that takes into account all the coupling among the nine correlation types. 
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B.3. Analytical Approximation 



The expressions for the coupling matrices greatly simplify when W pp > is diagonal in pixel 
space, Wp p i = 5pp>N bs,p- This limit is a good approximation to the WMAP data at high I, 
where noise is approximately uncorrelated (diagonal in pixel space). In this limit, one can 
evaluate the coupling matrices analytically. 



It is convenient to write the N f, s matrix as 



N QQ N QU 

obs,p obs,p 

N UQ N uu 

obs,p obs,p 



■^obs,p -^obs,p 
N UQ 
obs,p 



N QU 
obs,p 

~^obs,p ~^obs,p 



(B16) 



where 



obs,p 2 

atQQ _ atUU 
jy— obs,p obs,p 



(B17) 
(B18) 



obs,p 2 

One can show that under a rotation of basis by an angle 9, these quantities transform as 



NX 

obs,p 



Nt 

1 v o6s,p' 



Therefore, we expand them into spin harmonics as follows: 

^obs,p = n tm^lrn,p-, 
Im 

N : bs ,p ±iN oblp = J2^ nim ^ Yim - 



V 



(B19) 
(B20) 



(B21) 
(B22) 



Irn 



We obtain 
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LM 
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n 



LM 



EjLM 
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LM 

- [ + ,n LM + {-) L+l+l> 



AftLM 
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(B23) 



(B24) 
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4 n LM 



\L+l+V 



-4 n LM 



L I V 
-4 2 2 



(B31) 



C. Polarization Fisher and Covariance Matrix 



In this Appendix, we derive expressions for the Fisher and covariance matrices of the 
temperature and polarization power spec tra. Our derivation extends the derivation of the TT 
matrices given in iHinshaw et al.l (l2003bl ) to all combinations of polarization power spectra. 



Note that we do not use these results for evaluating the likelihood that is used in the 
cosmological analysis. At low multipoles, I < 23, we evaluate the likelihood of polarization 
data directly from the maps using the exact method described in Appendix [D] Why do we 
not use the Fisher or covariance matrix for the cosmological analysis, except for TT and TE 
spectra at I > 23? The reason is because the form of the likelihood function for the power 
spectra is not a Gaussian at low multipoles, and therefore the Fisher or covariance matrix, 
which only characterize the second-order moment of the power spectrum, is not sufficient 
to fully specify t he li kelihood function. This was pointed o ut after the first year release by 
Efstathioul (120041 ) and lSlosar et al.l (12004 ) and is discussed in lHinshaw et al.l (120061 ). As we do 
not know the precise form of the likelihood for the power spectra, we evaluate the likelihood 
of the temperature and polarization maps directly, which is a Gaussian, at low multipoles, 
/ < 23. For high multipoles, / > 23, the likelihood function may be approximated as a 
Gaussian and therefore we use a Gaussian likelihood with the Fisher or covariance matrices. 
While we do not use the EE or BB power spectra at I > 23, as they contain very little signal 
compared to noise, we do use the covariance matrix of the TE power spectrum at Z > 23 in 
the likelihood code, for which we adopt the analytical ansatz given in Equation IQ121 which 
was also used in the first- year analysis of the TE powe r spec trum (IKogut et al.ll2003l ). For 
the evaluation of the TT likelihood, see IHinshaw et al.l (120061 ). 



C.l. Fisher Matrix: Exact formula 



The Fisher matrix, i 7 ^/, is given by 



t^XY^'Y' 



In 



Be- 1 : 



dC. 



mi 9S xy 



(c- 1 ) 



dC, 



cm 
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where the covariance matrix C qq > consists of the covariance matrices of all the bilinear com- 
binations of T, Q, and U : 

/ r<TT /-iTQ n TU \ 

^r,W \ 



C qq i 



PP VP PP 

r QT r QQ r QU 

^pp' ^pp' ^pp' 

r<uT /Xjq r*uu 

\ ^pp' ^pp' ^pp' / 



(C2) 



and the covariance includes the signal and noise, C qq i = S qq > +N gq t. Here Sf Y is the angular 
(cross) power spectrum of the signal where X and Y denote T, E, or B. The inverse 
covariance matrix in harmonic space is then given by the harmonic transform of (C _1 ) £ 



iqq'- 



(C X )fm,Vm' - ^ ^(Xlm)(ap)(C *) ap ,a'p'Y ' (Y*l'm')(a'p') , (C3) 
ap,a'p' 

where T is given by the equations following (jB8J). 

Using these quantities, each term of the Fisher matrix (Eq. [CI] ) evaluates to 



t lV 


mm 1 


-uxx y 

Jim, I'm' J 
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(C4) 


r^XX,XY 
IV 


= EKc- 1 

mm' 




-yxv ] 

Jlm,l'm'\ ' 


(C5) 


r^XX,YY 
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mm 1 


-i\XY y 

Jim, I'm' J 


> 


(C6) 


TpXY,XY 


= EK^ 1 


>lm,l'm J ' 


" / J [(^ 1 )?m^'m'(^' 1 )hn^'m'] ' 


(C7) 



where 1^7. In general cases where S qq i or iV^/ (or both) are non-diagonal, one must 
calculate {C~ l ) qq i by directly inverting the covariance matrix given by equation ( 1C2I) . In 
reality, however, the matrix inversion requires rip operations and thus it become computa- 
tionally too expensive to evaluate for the full WMAP resolution. On the other hand, if one 
considers only large scale anisotropies at low I, then the matrix inversion can be done in a 
reasonable computational time. We use Eq. (1C4jl - (1C7jl for computing the Fisher matrices 
for Cf T , C[ E , Cf B , C EE , C EB , and Cf B , at low multipoles, I < 32. 



C.2. Fisher Matrix: Analytical Approximation 



The expressions for the Fisher matrices can be evaluated analytically when C qq > is diag- 
onal in pixel space. This limit is a good approximation to the WMAP data at high /, where 
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C qq t is dominated by noise and noise is approximately uncorrelated (diagonal in pixel space). 
In this limit, one obtains the following analytical formulae: 
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C.3. Covariance Matrix: Ansatz 

The inverse of the Fisher matrix gives the covariance matrix, S. While we use the 
map-based exact likelihood described in Appendix [D] for the cosmological analysis, it is still 
useful to have an approximate method to evaluate the likelihood of the data given theory 
and noise model from the power spectra. For this purpose, we use the following ansatz: 



^TETE 



(sF + <j s t ) (s™ + nff f t ) + (spy 



(2^+1) [/, 



TE 

sky eff 



(C12) 



-kTB tb 



(SJ T + n TT 



efft 



)(SF B + n°f fe ) 



(2£ + 1) [fF 



sky eff 
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(C13) 



^EEEE 



2(S EE + n? E e ) 



(2^ + 1) [fif 
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(C14) 



iBBBB 



2(Sf B + nff fe )'' 
(2£ + l\ \f BB 



IJ sky eff 



(C15) 



^EBEB 



(S EE + n EE 



efft 



)(Sr + n E f fe ) 



(C16) 

In these expressions n e ff£ denotes the effective noise as a function of £ and / s k ye // denotes 
the effective fraction of the sky observed. These are obtained from comparing the ansatz 
to the inverse of the Fisher matrices derived in the previous sections. We have found that 

fXY ^ 
J sky 



f sky f sky to a very good approximation. See also iKogut et al.l (120031 ) for the 



evaluation of T,J ETE and lHinshaw et al.l (120061 ) for the evaluation of S 



iTT TT 



D. Exact Likelihood Evaluation at Low Multipoles 

At low multipoles, I < 23, we evaluate the likelihood of the data for a given theoretical 
model exactly from the temperature and polarization maps. The standard likelihood is given 
by 

L(m\S)dm - avHrtiS + N)- 1 ™] dm 
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where rh is the data vector containing the temperature map, T, as well as the polarization 
maps, Q, and U, n p is the number of pixels of each map, and S and N are the signal and 
noise covariance matrix (3n p x 3n p ), respectively. As the temperature data are completely 
dominated by the signal at such low multipoles, noise in temperature may be ignored. This 
simplifies the form of likelihood as 



L(m\S)dm 



exp 


-\m\s P + Np) 


1 m 


dm ex P 1 


^-ifte/f) d f 




\Sp + Np\V* 




(2vr) n p 


\S T \ l l 2 (2tt)"p/2 



(D2) 



where St is the temperature signal matrix (n p x n p ), the new polarization data vector, 
m = (Q p , U p ), is given by 



u n = 



1 23 



Sf E 

1=2 1 m=-l 



23 rt 
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1=2 1 



/,ral+2^/m,p + -2^m,p)> 
lm{+2 1 lm,p ~ -2^ m ,p)' 



(D3) 
(D4) 



and Sp is the signal matrix for the new polarization vector with the size of 2n p x 2n p . As 

Ti m is totally signal dominated, the noise matrix for (Q, U) equals that for (Q, U), n p . 
To estimate Tim, we used the full-sky internal linear combination (LLC) temperature map 
jffinshaw et alJbnOfih . 

One can show that equation fIDll) and (ID2[) are mathematically equivalent when the 
temperature noise is ignored. The new form, equation (ID2I) . allows us to factorize the 
likelihood of temperature and polarization, with the information in their cross-correlation, 
Sf E , fully retained. We further rewrite the polarization part of the likelihood as 



L(m\S) 



exp 


-^{N^mYiN^SpNp 1 + Np 1 )- 1 ^? 1 ™) 


{N^ldm 




iN^SpNp 1 + Np 1 \ 1 / 2 
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(D5) 



This form is operationally more useful, as it contains only iVp 1 . iHinshaw et al.l (120061 ) 
describes the method to evaluate the temperature part of the likelihood. 



The effect of P06 mask is included in Np 1 . Suppose that the structure of Np 1 is given 



by 



Np 1 



A B 
B D 



(D6) 



where A is the noise matrix for unmasked pixels, D is for masked pixels, and B is for their 
correlations. We assign infinite noise to the masked pixels such that Np — > Np + X(I — 
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M), where M is the diagonal matrix whose elements are zero for masked pixels and unity 
otherwise. In the limit of A — > oo, the inverse of Np is given by 

! / A-&DB 
iVp ^ 

We have checked that this form of Np 1 yields the unbiased estimates of the signal matrix 
from simulated realizations of the WMAP data. When the masked pixels were simply ignored 
(i.e., B l DB = 0), on the other hand, the estimated signal matrix was found to be biased 
high. As the likelihood form is sensitive to the precise form of Np 1 , it is important to treat 
the mask in this way so that the estimated signal matrix from the data is unbiased. 

We mask the polarization maps as follows. We first mask the maps at the full resolution, 
inside = 512, and then degrade the masked maps using the weight that is diagonal in pixel 
space, Np 1 , to a lower resolution, n si d e = 16. (Note that while the weight is diagonal in 
pixel space, it contains noise covariance between Q p and U p . The spurious polarization term, 
S, is ignored in this process.) The degraded mask is redefined such that it takes on 1 when 
the lower resolution pixel contains more than half of the original full resolution pixels, and 
otherwise. We degrade these maps further to the resolution of n s ^ e = 8 using the full noise 
matrix, and also degrade the mask and the noise matrix. (The noise matrix has been masked 
using Eq |D7] .) We use the resulting maps and noise matrix in the likelihood function given 
in equation (1D5|) . 




E. An estimate of VLqw 

Tensor perturbations generated by inflation are stochastic in nature, so the gravity wave 
perturbation can be expanded in plane waves 

^3 [h^ V ,k)e±e- ik -* + h x ( V ,k)e* je - ik -*} , (El) 

where e"- is the polarization tensor, and a = +, x are the two polarizations in the transverse 
traceless (tt) gauge (in which hijj = h\ = 0; we also set h 00 = h j = 0). The stress-energy 
tensor for gravity waves is defined as 

^ = ^<WO, (E2) 

and in the tt gauge, we have 

T oo= 3^G<M«>. (E3) 



Thus, 
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The variance of the perturbations in the h fields can be written as 
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Writing 



h a (r),k) =h a (0,k)T(ri,k), 
where T is the transfer function, we have 
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From the definition of the primordial tensor power spectrum, 
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Remembering that f2 = p x (8ttG/3Hq), we obtain 
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Therefore, 



n 



GW 



dink 



12H% 



(E14) 



The transfer function T and its time derivative t can be calculated easily by numerically 
integrating the evolution equation for the polarization states, which, neglecting the neutrino 
anisotropic stress, is given by 



K + 2 - 
a 



K + k 2 h a = o, 



(E15) 



where prime denotes derivatives with respect to conformal time r), related to the time deriva- 
tive by dr] = dt/a(rj). This expression may be numerically integrated. In the following, 
however, we derive an analytic estimate relating a given limit on the tensor-to-scalar ratio, 
r, and the measured amplitude of the primordial scalar power spectrum, A, to a limit on the 
current energy density in primordial gravitational radiation. 

There are several approaches taken in the literature to derive analytic expressions for the 
tensor transfer function, though these results are obtained in almost all cases for a universe 
containing only matter and radiation. These include using (1) an instantaneous transition 
from radiation to matter domination jAbbott fc Harari 1986 ; Ng fc Speliotopoulos 1995 



Grishchukll200ll ; |Pritchard fc Kamionkowskill2005l . e.g., ) (2) a "transfer functio n" to account 



1993 



Wang 



1996 



Turner 



1997 



Pritchard fc Kamionkowski 



for the smooth transition from radiation domination to matte r domination (ITurner et al. 

e.g., ), and (3) WKB methods JNg fc Speliotopoulosl fl9951 : 
20051 . e.g., ). In the foll owing derivation, we will apply the 



sudden transition approximation to a ACDM universe (jZhang et al.ll2005l . see also), which is 
a good approximation for gravitational waves with wavelengths much longer than the time 
taken for the transition to happen. 

In a universe which undergoes a set of piecewise instantaneous transitions in the scale- 
factor, given by a{rf) oc r]~ u , the solution to eq. IE15I is given by 

(k v y +1 [Cj u (k V )+D yu (krj)}, 



h(rj,k) 



(E16) 



where j u and y v are spherical Bessel functions of order v of the first and second kinds, 
respectively. Here, v = — 1 for radiation-domination (RD, 77 < T] eq i), v = —2 for matter- 
domination (MD, i] eq i < r] < r] eq 2), and v = +1 for A-domination (LD, r\ > 77^2) • Veqi 
is the conformal time at radiation-matter equality, with a scale-factor corresponding to 
(Q r /Q m ), and rj eq2 is the conformal time at matter-A equality, with a scale-factor 



L eql 



corresponding to a, 



eq2 



1 1/3 



For a concordance cosmology with {Q r , fl m , Q\, h} 



{4.18 x 10- 5 //i 2 , 0.3, 0.7, 0.72}, r] eql = 103 Mpc" 1 and r} eg2 
Mpc^ 1 respectively in the instantaneous approximation). 



12270 Mpc" 1 (115 and 12030 
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To obtain the coefficients C and D, we require h and h' to be continuous at each of the 
transitions, r) eql and i] eq2 . Thus, denoting x = ki] and making use of special properties of 
spherical Bessel functions, we obtain the transfer function and its derivative at present: 

T(x) = x 2 [C j 1 (x) + Dy 1 (x)}, (E17) 
f(k,x) = kx 2 [C j (x) + Dy (x)}. (E18) 

The coefficients are given by 

C = [2Ax\ + 35(1 + x 2 2 ) + 3 cos(2:r 2 ) (B + 2Ax 2 - Bxf) + 3 sin(2:r 2 ) {2Bx 2 + A{x\ - 1))] 

(E19) 

D = -L [2Bx\ - 3,4(1 + ^2) + 3 cos(2a; 2 ) (A - 2Bx 2 - Ax 2 2 ) + 3 sin(2:r 2 ) (B + 2Ax 2 - Bxj)] 

(E20) 

3xi - XiC0s(2xi) +2sin(2xi) 

^ = 2^ (E21) 

B = 2-2x\-2 cos(2a: 1 ) - x x anggQ 22 

where xi = fc?7 eg i and x 2 = fc?7 eg2 . 

Further, we have the following definitions: 

\ «t(feo) 



A^(fc) = A K k »)[f) (E23) 

' ■ »■ ° 

where 

A^(A; ) ~ 2.95 x KT 9 A(k ). (E25) 
To eliminate n t , we use the inflationary single-field consistency relation, n t = — r/8. 

Combining these equations, and evaluating them at the present conformal time i] (with 
a = 1) for modes within our current horizon, we are left with 

1 f°° dk ( k \ ~ r/8 2 

Q gw = T77772 / ~T rA n( k o) t k2 ( k Vo) 4 [C(k, rj eql , rj eq2 )j (kri ) + D(k, rj eql , r] eq2 )y (kr] )} , 

(E26) 

where k and 770 are to be evaluated in units of 1/Mpc and k = 0.002 Mpc -1 . We can now 
change to the dimensionless variable x = krjo and obtain 

n G w - 2,95 X ^ r dx x "~ r/8 [^o(x) + A/o(*)] 2 , (E27) 

izn r] j 2n 
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where xq = k rj . We also have the result 



GW 



dink 



(k, r)o) = 2.21 x 10- 3 (rA) - 



-r/8 



k 

k N - r/8 



T{k,r]o) 



-i 2 



(E28) 



2.21 x KT 3 (rA) ( - ) {kx 2 [C Jo (x) + D y {x)]y . (E29) 



Now 



H r] 



(E30) 



y/Q r + Vt m a + ^A a4 ' 

and -ffo^o — 3.25 for the concordance ACDM model. Taking the concordance model and 
fcn = 0.002 M pc" 1 , for given upper limits on r and A, the upper limit on Vt GW is given by 
JPeirisll2003h . 

tt GW < 2.33xl0" n (r A) (27.05) r/8 [0.1278 - 0.0835 (logr) - 0.0671 (logr) 2 - 0.0248 (logr) 3 ] 

(E31) 

where the logarithm is taken in base ten. 



REFERENCES 

Abbott, L. F. & Harari, D. D. 1986, Nucl. Phys., B264, 487 

Baccigalupi, C, Burigana, C, Perrotta, F., De Zotti, G., La Porta, L., Maino, D., Maris, 
M., & Paladini, R. 2001, A&A, 372, 8 

Banday, A. J. & Wolfendale, S. 1991, MNRAS, 248, 705 

Barkats, D. et al. 2005, ApJ, 619, L127 

Barnes, C, et al. 2002, ApJS, 143, 567 

-. 2003, ApJS, 148, 51 

Basko, M. M. & Polnarev, A. G. 1980, MNRAS, 191, 207 

Beck, R. 2001, Space Science Reviews, 99, 243, kluwer Academic Publishers 

Beck, R. 2006, in Polarisation 2005, ed. F. B. . M. Miville-Deschenes (EAS Publication 
Series) 

Beck, R., Brandenburg, A., Moss, D., Shukurov, A., & Sokoloff, D. 1996, Annu. Rev. Asron. 
Astrophys., 34, 155 



Bennett, C. L., et al. 2003a, ApJS, 148, 97 

-. 2003b, ApJ, 583, 1 

Benoit, A., et al. 2004, A&A, 424, 571 

Berdyugin, A., Piirola, V., & Teerikorpi, P. 2004, A&A, 424, 873 
Berdyugin, A. & Teerikorpi, P. 2001, a, 368, 635 
-. 2002, A&A, 384, 1050 

Berdyugin, A., Teerikorpi, P., Haikala, L., Hanski, M., Knude, J., & Markkanen, T. 2001, 
A&A, 372, 276 

Berkhuijsen, E. M., Haslam, C. G. T., & Salter, C. J. 1971, A&A, 14, 252 

Bernardi, G., Carretti, E., Fabbri, R., Sbarra, C, Poppi, S., & Cortiglioni, S. 2003, MNRAS, 
344, 347 

Blum, E. J. 1959, Annales D' Astrophysique, 22, 140 
Bond, J. R. & Efstathiou, G. 1984, ApJ, 285, L45 
Bond, J. R. & Efstathiou, G. 1987, MNRAS, 226, 655 

Born, M. & Wolf, E. 1980, Principles of Optics, sixth edn. (Pergamon Press) 

Boyle, L., Steinhardt, P., & Turok, N. 2005, Phys. Rev. Lett., Phys. Rev. Lett., submitted, 
( |astro-ph/0507455l ) 

Brouw, W. N. & Spoelstra, T. A. T. 1976, Astron. Astrophys. Suppl. Ser., 26, 129 

Bunn, E. F., Zaldarriaga, M., Tegmark, M., & de Oliveira-Costa, A. 2003, Phys. Rev. D, 67, 
23501 

Burn, B. J. 1966, MNRAS, 133, 67 

Caderni, N., Fabbri, R., Melchiorri, B., Melchiorri, F., & Natale, V. 1978, Phys. Rev. D, 17, 
1908 

Cioffi, D. F. & Jones, T. W. 1980, AJ, 85, 368 
Cortiglioni, S. & Spoelstra, T. A. T. 1995, A&A, 302, 1 

Coulson, D., Crittenden, R. G., & Turok, N. G. 1994, Phys. Rev. Lett., 73, 2390 



-61 - 

Crittenden, R., Davis, R. L., & Steinhardt, P. J. 1993, ApJ, 417, L13 
Crittenden, R. C, Coulson, D., & Turok, N. G. 1995, Phys. Rev. D, 52, 5402 
Davis, L. J. & Greenstein, J. L. 1951, ApJ, 114, 206 

de Oliveira-Costa, A., Kogut, A., Devlin, M. J., Netterfield, C. B., Page, L. A., & Wollack, 
E. J. 1997, ApJ, 482, L17+ 

de Oliveira-Costa, A., Tegmark, M., Page, L., & Boughn, S. 1998, ApJ, 509, L9 

de Oliveira-Costa, A. et al. 1999, ApJ, 527, L9 

-. 2003, Phys. Rev. D, 68 

-. 2004, Astrophys. J., 606, L89 

Dodelson, S. 2003, Modern Cosmology (Academic Press) 
Draine, B. T. & Lazarian, A. 1998, ApJ, 494, L19 
-. 1999, ApJ, 512, 740 

Drimmel, R. & Spergel, D. N. 2001, ApJ, 556, 181 

Dumke, M., Krause, M., Wielebinski, R., & Klein, U. 1995, Astro. & Astro., 302, 691 
Duncan, A., Reich, P., Reich, W., & Furst, E. 1999, A&A, 350, 447 
Duncan, A. R., Haynes, R. F., Jones, K. L., & Stewart, R. T. 1995, MNRAS, 277, 36 
Efstathiou, G. 2004, MNRAS, 348, 885 

Eisenhauer, F., Schoedel, R., Genzel, R., Ott, T., Tecza, M., Abuter, R., Eckart, A., & 
Alexander, T. 2003, ApJ, 597, 121 

Enomoto, R., et al. 2002, Nature, 416, 823 

EnBlin, T. A., Waelkens, A., Vogt, C, & Schekochihin, A. A. 2006, Astronomische 
Nachrichten, 327, 626 

Erickson, W. C. 1957, ApJ, 126, 480 

Faris, J. 1967, J. Res. Nat. Bur. Stand., Engr. & Instr., 71C, 153 
Fernandez- Cerezo, S., et al. 2006, ArXiv Astrophysics e-prints 



-62 - 

Finkbeiner, D. P. 2003, ApJS, 146, 407, accepted ( |astro-ph/030"l558| 
-. 2004, ApJ, 614, 186 

Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 

Finkbeiner, D. P., Langston, G. I., & Minter, A. H. 2004, ApJ, 617, 350 

Finkbeiner, D. P., Schlegel, D. J., Frank, C, & Heiles, C. 2002, ApJ, 566, 898 

Flett, A. M. & Henderson, C. 1979, MNRAS, 189, 867 

Flett, A. M. & Murray, A. G. 1991, MNRAS, 249, 4P 

Fosalba, P., Lazarian, A., Prunet, S., & Tauber, J. A. 2002, ApJ, 564, 762 

Frewin, R. A., Polnarev, A. G., & Coles, P. 1994, MNRAS, 266, L21+ 

Giardino, G., Banday, A. J., Gorski, K. M., Bennett, K., Jonas, J. L., & Tauber, J. 2002, 
A&A, 387, 82 

Gorski, K. M., Hivon, E., & Wandelt, B. D. 1998, in Evolution of Large-Scale Structure: 
From Recombination to Garching 

Greaves, J. S., et al. 2003, MNRAS, 340, 353 

Grishchuk, L. P. 2001, Lect. Notes Phys., 562, 167 

Han, J.-L. 2006, Chinese Journal of Astronony and Astrophysics, submitted 
dasfao-ph/0603512"l ) 



Han, J.-L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van Straten, W. 2006, ApJ, 
accepted ( |astro-ph/ 0601357) 

Harari, D. D. & Zaldarriaga, M. 1993, Physics Letters B, 319, 96 

Haslam, C. G. T., Stoffel, H., Salter, C. J., & Wilson, W. E. 1982, A&AS, 47, 1 

Hedman, M. M., Barkats, D., Gundersen, J. O., McMahon, J. J., Staggs, S. T., & Winstein, 
B. 2002, ApJ, 573, L73 

Heiles, C. 2000, AJ, 119, 923 

Hinshaw, G., et al. 2003a, ApJS, 148, 63 

-. 2003b, ApJS, 148, 135 



-63- 

Hinshaw, G. et al. 2006, ApJ, submitted 

Hivon, E., Gorski, K. M., Netterfield, C. B., Crill, B. P., Prunet, S., & Hansen, F. 2002, 
ApJ, 567, 2 

Hobbs, R. W., Maran, S. P., & Brown, L. W. 1978, ApJ, 223, 373 
Holder, G. P., Haiman, Z., Kaplinghat, M., & Knox, L. 2003, ApJ, 595, 13 
Hu, W., Hedman, M. M., & Zaldarriaga, M. 2003, Phys. Rev. D, 67, 43004 
Hu, W. & White, M. 1997, Phys. Rev. D, 56, 596 

Hummel, E., Dahlem, M., van der Hulst, J. M., & Sukumar, S. 1991, A&A, 246, 10 

Jarosik, N., et al. 2003a, ApJS, 145, 413 

-. 2003b, ApJS, 148, 29 

Jarosik, N. et al. 2006, ApJ, submitted 

Johnston, K. J. & Hobbs, R. W. 1969, ApJ, 158, 145 

Jones, R. C. 1941, J. Opt. Soc. Am., 31, 488 

Jones, T. J., Klebe, D., & Dickey, J. M. 1992, ApJ, 389, 602 

Junkes, N., Haynes, R. F., Harnett, J. I., & Jauncey, D. L. 1993, A&A, 269, 29 

Kaiser, N. 1983, MNRAS, 202, 1169 

Kamionkowski, M., Kosowsky, A., & Stebbins, A. 1997, Phys. Rev. D, 55, 7368 

Keating, B. G., O'Dell, C. W., de Oliveira-Costa, A., Klawikowski, S., Stebor, N., Piccirillo, 
L., Tegmark, M., & Timbie, P. T. 2001, ApJ, 560, LI 

Keating, B. G., Timbie, P. T., Polnarev, A., & Steinberger, J. 1998, ApJ, 495, 580 

Khoury, J., Ovrut, B. A., Seiberg, N., Steinhardt, P. J., & Turok, N. 2002, Phys. Rev., D65, 
086007 

Kogut, A., Banday, A. J., Bennett, C. L., Gorski, K. M., Hinshaw, G., & Reach, W. T. 
1996, ApJ, 460, 1 

Kogut, A., et al. 2003, ApJS, 148, 161 



-64- 



Kosowsky, A. 1996, Ann. Phys., 246, 49 

Kovac, J. M., Leitch, E. M., Pryke, C, Carlstrom, J. E., Halverson, N. W., k Holzapfel, 
W. L. 2002, Nature, 420, 772 

Krauss, L. M. k White, M. 1992, Physical Review Letters, 69, 869 

Lagache, G. 2003, A&A, 405, 813 

Lawson, K. D., Mayer, C. J., Osborne, J. L., k Parkinson, M. L. 1987, MNRAS, 225, 307 
Lazarian, A. k Draine, B. T. 2000, ApJ, 536, L15 

Leitch, E. M., Kovac, J. M., Halverson, N. W., Carlstrom, J. E., Pryke, C, k Smith, 
M. W. E. 2005, ApJ, 624, 10 

Leitch, E. M., Readhead, A. C. S., Pearson, T. J., k Myers, S. T. 1997, ApJ, 486, L23 

Leitch, E. M., et al. 2002, Nature, 420, 763 

Lewis, A., Challinor, A., k Lasenby, A. 2000, ApJ, 538, 473 

Lewis, A., Challinor, A., k Turok, N. 2002, Phys. Rev. D, 65, 023505 

Liddle, A. R. k Lyth, D. H. 2000, Cosmological inflation and large-scale structure (Cam- 
bridge University Press) 

Limon, M., et al. 2003, Wilkinson Microwave Anisotropy Probe (WMAP): Explanatory 
Supplement, http : //lambda, gsfc .nasa. gov/data/map/doc/MAP_supplement .pdf 

Linde, A. D. 1983, Phys. Lett., B129, 177 

Lubin, P., Melese, P., k Smoot, G. 1983, ApJ, 273, L51 

Lubin, P. M. k Smoot, G. F. 1979, Physical Review Letters, 42, 129 

-. 1981, ApJ, 245, 1 

Matveenko, L. I. k Conklin, E. K. 1973, Soviet Astronomy, 16, 726 
Mayer, C. H. k Hollinger, J. P. 1968, ApJ, 151, 53 
Mitner, A. H. k Spangler, S. R. 1996, ApJ, 458, 194 

Montgomery, C. G., Dicke, R. H., k Purcell, E. M. 1948, MIT Radiation Laboratory Series, 
Vol. 8, Principles of Microwave Circuits (New York: McGraw-Hill) 



-65 - 



Montgomery, J. W., Epstein, E. E., Oliver, J. P., Dworetsky, M. M., & Fogarty, W. G. 1971, 
ApJ, 167, 77 

Montroy, T. E., et al. 2005, ApJ 

Mukhanov, V. 2005, Physical Foundations of Cosmology (Cambridge University Press) 
Mukhanov, V. & Vikman, A. 2005, LMU-ASC 78/05 

Mukherjee, P., Coble, K., Dragovan, M., Ganga, K., Kovac, J., Ratra, B., & Souradeep, T. 
2003, ApJ, 592, 692 

Nanos, G. P. 1979, ApJ, 232, 341 

Netterfield, C. B., Devlin, M. J., Jarosik, N., Page, L., & Wollack, E. J. 1997, ApJ, 474, 47 

Ng, K. L. & Ng, K.-W. 1995, Phys. Rev., D51, 364 

Ng, K.-W. & Speliotopoulos, A. D. 1995, Phys. Rev., D52, 2112 

Page, L., et al. 2003a, ApJS, 148, 39 

-. 2003b, ApJ, 585, 566 

Peebles, P. J. E. 1968, ApJ, 153, 1 

— . 1993, Principles of Physical Cosmology (Princeton, NJ: Princeton University Press) 

Peiris, H. 2003, Ph.D. thesis, Princeton University 

Peiris, H. V., et al. 2003, ApJS, 148, 213 

Penzias, A. A. & Wilson, R. W. 1965, ApJ, 142, 419 

Piacentini, F., et al. 2005, ArXiv Astrophysics e-prints 

Polnarev, A. G. 1985, AZh, 62, 1041 

Ponthieu, N. et al. 2005, Astron. & Astro,, 607, 655 

Pritchard, J. R. & Kamionkowski, M. 2005, Annals Phys., 318, 2 

Prunet, S. & Lazarian, A. 1999, in Microwave Foregrounds. Sloan Summit, Institute for 
Advanced Study, Princeton, New Jersey 14-15 November 1998. Eds: A. Rde Oliveira- 
Costa and M. Tegmark., Vol. 181 (Astronomical Society of the Pacific), 113 



-66- 



Readhead, A. C. S., et al. 2004, ApJ, 609, 498 
Readhead, A. C. S., et al. 2004, Science, 306, 761 
Rees, M. J. 1968, ApJ, 153, LI 
Reich, R & Reich, W. 1988, A&AS, 74, 7 

Reich, W. 2006, in Cosmic Polarization, Ed: Roberto Fabbri (Research Signpost) 

Reid, M. & Brunthaler, A. 2005, in Future Directions in High Resolution Astronomy: The 
10th Anniversary of the VLBA, ASP Conference Proceedings, Eds: J. Romney and 
M. Reid., Vol. 340 (Astronomical Society of the Pacific), 253 

Rybicki, G. B. & Lightman, A. 1979, Radiative Processes in Astrophysics (Wiley & Sons: 
New York) 

Sault, R. J., Hamaker, J. P., & Bregman, J. D. 1996, A&AS, 117, 149 

Seljak, U. 1997, ApJ, 482, 6 

Seljak, U. & Zaldarriaga, M. 1996, ApJ, 469, 437 

Sievers, J. L., et al. 2005, ArXiv Astrophysics e-prints 

Sironi, C, Boella, C, Bonelli, C, Brunetti, L., Cavaliere, F., Gervasi, M., Giardino, G., & 
Passerini, A. 1997, New Astronomy, 3, 1 

Slosar, A., Seljak, U., & Makarov, A. 2004, Phys. Rev., D69, 123003 

Sofue, Y., Fujimoto, M., & Wielebinski, R. 1986, ARA&A, 24, 459 

Spergel, D. N. & Zaldarriaga, M. 1997, Phys. Rev. Lett., 79, 2180 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Spergel, D. N. et al. 2006, ApJ, submitted 

Starobinsky, A. A. 1979, ZhETF Pis ma Redaktsiiu, 30, 719 

Steinhardt, P. J. & Turok, N. 2002, Science, 296, 1436 

Sukumar, S. & Allen, R. 1991, ApJ, 382, 100 

Taylor, J. H. & Cordes, J. M. 1993, ApJ, 411, 674 



-67- 

Tegmark, M., Eisenstein, D. J., Hu, W., & de Oliveira-Costa, A. 2000, ApJ, 530, 133, 
asto>ps/9905257| 

Tinbergen, J. 1996, Astronomical polarimetry (Cambridge University Press) 
Torbet, E. et al. 1999, ApJ, 521, L79 

Tucci, M., Carretti, E., Cecchini, S., Nicastro, L., Fabbri, R., Gaensler, B. M., Dickey, J. M., 
& McClure-Griffiths, N. M. 2002, ApJ, 579, 607 

Turner, M. S. 1997, Phys. Rev. D, 55, 435 

Turner, M. S., White, M. J., & Lidsey, J. E. 1993, Phys. Rev., D48, 4613 
Turok, N. 1996, ApJ, 473, L5 

Uyaniker, B., Fuerst, E., Reich, W., Reich, P., & Wielebinski, R. 1999, Astron. Astrophys. 
Suppl. Ser., 138, 31 

Verde, L., Peiris, H., & Jimenez, R. 2006, JCAP, 019 

Verde, L., et al. 2003, ApJS, 148, 195 

Wang, Y. 1996, Phys. Rev., D53, 639 

Watson, R. A., Rebolo, R., Rubino-Martm, J. A., Hildebrandt, S., Gutierrez, C. M., 
Fernandez-Cerezo, S., Hoyland, R. J., & Battistelli, E. S. 2005, ApJ, 624, L89 

Weiss, R. 1984, ARA&A, 12, 90 

Wielebinski, R. 2005, in Cosmic Magnetic Fields, ed. R. W. . R. Reck, Vol. 664 (LSP, Berlin: 
Springer), 89 

Wollack, E. B., Jarosik, N., Netterfield, C. B., Page, L., & Wilkinson, D. 1993, ApJ, 419, 
L49 

Wolleben, M., Landecker, T. L., Reich, W., & Wielebinski, R. 2005, ArXiv Astrophysics 
e-prints 

Wright, E. L., Hinshaw, G., & Bennett, C. L. 1996, ApJ, 458, L53 
Wright, M. C. H. & Forster, J. R. 1980, ApJ, 239, 873 
Zaldarriaga, M. & Harari, D. D. 1995, Phys. Rev., D52, 3276 



Zaldarriaga, M. & Seljak, U. 1997, Phys. Rev. D, 55, 1830 

Zeldovich, Y. B., Kurt, V. G., & Syunyaev, R. A. 1969, Journal of Experimental and Theo- 
retical Physics, 28, 146 

Zhang, Y., Yuan, Y., Zhao, W., & Chen, Y.-T. 2005, Class. Quant. Grav., 22, 1383 



This preprint was prepared with the AAS IATgX macros v5.2. 



-69- 



Table 1. Polarization of Tau A 



Measurement 


v [GHz] 


i [Jy] 


Q [Jy] 


U [Jy] 


P/I [%] 


jpa [deg] 


WMAP 


22.5 (K) 


352 ± 11 


-24.7 ±0.8 


1.3 ±0.9 


7.0 ±0.3 


-88? (150°) 


WMAP 


32.8 (Ka) 


322 ±6 


-22.2 ±2.0 


1.9 ± 1.1 


6.9 ±0.3 


-87? (151°) 


WMAP 


40.4 (Q) 


299 ±6 


-19.6 ±2.6 


0.5 ±2.4 


6.6 ±0.9 


-89? (149°) 


WMAP 


60.2 (V) 


265 ±7 


-18.5 ±2.7 


-1.9 ±6.2 


7.0 ± 1.1 


-93? (145°) 


WMAP 


92.9 (W) 


229 ±11 


-17.5 ±4.4 


-1.3 ±7.2 


7.6 ±2.0 


-92? (146°) 


Maver & Hollineer (1968) 


19 








6.6 [15.5] 


(140° ± 10) 


Wright & Forster (1980) 


23 








9 


(152°) 


Johnston & Hobbs (19691 


31 








8.1 [17] 


(158°) 


Flett & Henderson (19791 


33 








[16] 


([154?8±2]) 


Matveenko & Conklin (19731 


8G 








([23 ±3]) 




Montgomery et al. (19711 


88 








13 


(152°) 


Hobbs et al. (19781 


99 








[11.9 ±0.9] 


([123°]) 


Flett & Murray (19911 


273 








[27 ± 1] 


([146° ±2]) 


Greaves et al. (20031 


363 








25 ±5 


(150° ±6) 



Note. - The fluxes are integrated over pixels within a radius that includes 99% of the beam solid 
angle, r99 = [2?525, 1?645, 1?517, 1?141,0?946] degrees in K through W bands. The errors are la estimates 
calculated as a quadrature sum of statistical error, error due to background uncertainty, confusion error, 
0.5% calibration error, and an additional 1% error since the aperture radius does not include all of the beam 
solid angle. Confusion error was calculated as the maximum difference in derived flux when the aperture 
radius and annulus radius are both decreased by 20% or increased by 20%. Confusion error is usually the 
largest contribution to the total error. The frequencies are band center frequencies for Tau A's antenna 
temperature spectral index, (3 = —2.3. The two numbers for jpa correspond to Galactic and equatorial 
(in parentheses) coordinates. Non- WMAP measurements are generally done with arcminute resolution and 
therefore have different average and peak (in square brackets) fractional polarization. Their polarization 
directions are all in equatorial coordinates. 
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Table 2. Temperatures in the Galactic Center Region 



Band 


I [mK] 


Q [mK] 


U [mK] 


K 


33 


0.69 


-0.25 


Ka 


14 


0.21 


-0.086 


Q 


8.7 


0.10 


-0.041 


V 


4.0 


0.037 


-0.01 < U < 0.01 


w 


3.6 


0.043 


-0.01 < U < 0.01 



Note. - The table gives the average values for the temperature and Q and U Stokes 
parameters in a 5b = 2° by SI = 10° region centered on (l,b) = (0,0). The values are in 
thermodynamic units relative to the CMB. To convert to antenna temperature, divide by 
1.014, 1.029, 1.044, 1.100, 1.251 in K through W bands respectively. 
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Table 3. Fit Coefficients to Foreground Templates 



Band 






OLd,v 




Ka 


0.3103 


-3.22 


0.0148 


1.54 


Q 


0.1691 


-3.12 


0.0154 


1.89 


V 


0.0610 


-2.94 


0.0343 


1.92 


w 


0.0358 


-2.51 


0.0891 




Ka 


0.2973 


-3.33 


0.0148 


1.54 


Q 


0.1492 


-3.33 


0.0154 


1.89 


V 


0.0414 


-3.33 


0.0343 


1.92 


w 


0.0112 


-3.33 


0.0891 





Note. - The top of the Table gives the coefficients for a direct fit to the polarization 
maps. The a are dimensionless and produce model maps in thermodynamic units. The 
spectral indices (3 refer to antenna temperature. The bottom half of the Table gives the 
same numbers for when the synchrotron fit is constrained to follow a power law. The fits 
were evaluated outside the processing mask. 
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Table 4. 


Comparison of x 2 


Between Pre-cleaned and cleaned Maps 


Band 


X 2 /v Pre-cleaned X 2 / u Cleaned v 


Ax 2 


Ka 


10.65 


1.20 


6144 


58061 


Q 


3.91 


1.09 


6144 


17326 


V 


1.36 


1.19 


6144 


1045 


w 


1.38 


1.58 


6144 


-1229 


Ka 


2.142 


1.096 


4534 


4743 


Q 


1.289 


1.018 


4534 


1229 


V 


1.048 


1.016 


4534 


145 


w 


1.061 


1.050 


4534 


50 



Note. — The top half of the table compares x 2 l v f°r the full-sky pre-cleaned map to x 2 l v 
for full-sky cleaned map. The bottom half makes a similar comparison for the region outside 
the P06 mask. 
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lable 5 


. WMAP EE i=2 and BB^ = 


=5 Values for £(£ + l)C e /2n 






Cross 


Bin 


feff 




yl-y2 EE^ =2 


BB^ =5 


yl-y2 BB e=5 


j^^C'leaned 








GHz 


(^K) 2 


(^K) 2 


(^K) 2 


(/.K) 2 


(^K) 2 


(/"K) 2 


KK 


1 


22.8 


Qnfi a 4- n 1 o 
oUd.d ± 




oo.l ± U.lo 








KKa 


2 


27.4 


93.3 ± 0.07 


0.0 ±0.22 


14.9 ±0.10 


0.6 ±0.30 


0.8 ±0.10 


0.6 ±0.14 


KQ 


2 


30.5 


53.6 ± 0.09 


-1.6 ±0.27 


8.5 ±0.11 


-1.3 ±0.32 


o r\ i n 1 n 

3.0 ± 0.10 


0.2 ± 0.11 


KV 


3 


37.2 


21.8 ±0.10 


-0.7 ±0.29 


2.0 ±0.13 


-0.6 ±0.38 


l.6±0.l0 


-0.7 ±0.13 


KW 




46.2 


10.4 ±0.13 


-3.8 ±0.4 


0.1 ±0.17 


-0.3 ±0.52 


-7.4 ±0.14 


-1.9 ±0.18 


KaKa 


2 


33.0 


30.5 ±0.13 




4.8 ±0.17 




0.7 ±0.26 


-0.1 ±0.35 


KaQ 


3 


36.6 


17.2 ± 0.09 


-0.0 ±0.27 


2.7 ±0.11 


-0.7 ±0.32 


0.6 ±0.15 


-0.1 ±0.18 


KaV 


4 


44.8 


8.2 ±0.10 


0.2 ±0.30 


0.7 ±0.12 


0.2 ±0.37 


0.1 ±0.15 


-0.2 ±0.19 


KaW 


4 


55.5 


5.9 ±0.14 


0.6 ±0.41 


0.6 ±0.17 


0.0 ±0.51 


0.4 ±0.20 


-0.1 ±0.25 


QQ 


4 


40.7 


9.6 ±0.17 


-0.1 ±0.67 


1.8 ±0.17 


0.3 ±0.68 


0.3 ±0.23 


0.0 ±0.24 


QV 


4 


49.7 


4.5 ±0.12 


-0.1 ±0.37 


0.6 ±0.13 


0.9 ± 0.40 


-0.1 ±0.15 


0.0 ±0.16 


QW 


4 


61.7 


3.3 ±0.17 


0.2 ±0.5 


0.7 ±0.18 


-0.1 ±0.55 


0.1 ±0.20 


0.2 ±0.21 


vv 


4 


60.8 


2.4 ±0.21 


-0.5 ±0.81 


0.2 ±0.21 


-0.2 ±0.65 


0.5 ±0.19 


0.2 ±0.23 


vw 


4 


75.4 


2.3 ±0.18 


1.0 ±0.55 


0.2 ±0.21 


-0.2 ±0.65 


0.5 ±0.19 


0.2 ±0.23 


WW 


4 


93.5 


2.2 ±0.37 


1.5 ± 1.27 


-0.4 ±0.44 


-0.3 ± 1.48 


0.3 ±0.38 


-0.7 ±0.45 



Note. — For v > 40 GHz, the largest foreground signals are at t = 2 of EE and £ = 5 of BB. 
This table shows the "raw" and "cleaned" values. The column labeled "bin" indicates which cross 
spectra are coadded into frequency bins. Because K band is used as a foreground template, there 
are no foreground corrected values. Also, as there are only single K and Ka band polarization 
channels, it is not possible to form cross spectra of year one minus year two. The yl-y2 notation 
refers to year one minus year two. KW is not used in any of the averages over frequency. 
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Table 6. x 2 /" for TE, TE (2003), TB, EE, BB, and EB 





r4 a 


r gb,f 


i-9 


r9 




I = 2 - 16 


I = 17-100 


I = 17-500 


1 = 17- 800 




(u = 15 dof) 


(v = 84 dof) 


colhead(z/ = 484 dof) 


{y = 784 dof) 


TE c ' e,g,h 


0.31 (0.99) f 


1.01 (0.46) 


1.20(0.01) 


1.08(0.06) 


TE (2003) d 


1.88 (0.03) 


1.18 (0.25) 


2.06 (0) 




XB c,g 


0.57 (0.90) 


0.72 (0.97) 


0.97 (0.70) 


0.97 (0.74) 


EE C 


1.34 (0.17) 


1.06 (0.33) 


0.98 (0.59) 


0.96 (0.76) 


BB C 


0.72 (0.77) 


1.28 (0.04) 


0.96 (0.73) 


0.95 (0.81) 


EB C 


0.41 (0.98) 


1.21 (0.09) 


1.03 (0.34) 


0.96 (0.76) 



Note. — x l v is computed for the null model {Cf = 0). 

a r4 HEALPix maps are used for i < 32. We limit this to I < 17 to avoid pixel window 
effects. 

b r9 HEALPix maps are used for 16 < £ < 800. 

c For all results a model of the foreground emission has been removed. 



d TE (2003) corresponds to lKogut et all (120031 ) . 
e The numbers in parentheses are the PTEs. 

f For £ > 16, we use the binned diagonal elements of the covariance matrices in Ap- 
pendix IC.31 

g For TE and TB, the E and B are comprised of a combination of Q and V bands and the 
T is from V and W bands. 

h The TE signal is in this I range and so the PTE should be low. 
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Table 7. x 2 1 v f° r r 4 Yearly Difference Null Maps 





yi-y2 


y2-y3 


yi - y3 




/ = 2 - 16 (i/ = 15 dof) 


I = 2 - 16 (u = 15 dof) 


l = 2 - 16 (u = 15 dof) 


TE 


1.70 (0.04) 


1.05 (0.40) 


1.87 (0.02) 


TB 


1.95 (0.02) 


1.20 (0.26) 


1.08 (0.37) 


EE 


1.55 (0.08) 


0.89 (0.58) 


0.55 (0.91) 


BB 


0.56 (0.90) 


1.50 (0.09) 


0.76 (0.72) 


EB 


0.62 (0.86) 


1.04 (0.41) 


0.84 (0.63) 



Note. — x 2 l v is computed for the null model, Cf x = 0. 
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Table 8. Binned Data for B EE /£ for £ > 20 





30 < I < 50 


51 < e < 150 


151 < t < 250 


251 < I < 350 


351 < I < 450 


451 < t < 650 


651 < t < 1023 


QV 


0.010 ± 0.007 


0.011 ± 0.005 


-0.001 ± 0.012 


-0.003 ± 0.026 


-0.014 ± 0.058 


0.16 ± 0.12 


-0.73 ± 0.66 


vw 


0.013 ± 0.011 


0.004 ± 0.004 


0.017 ± 0.009 


0.027 ± 0.018 


0.031 ± 0.037 


0.095 ± 0.065 


0.13 ± 0.22 


QVW 


0.013 ± 0.006 


0.004 ± 0.004 


0.017 ± 0.009 


0.027 ± 0.018 


0.031 ± 0.037 


0.095 ± 0.065 


0.13 ± 0.22 


KaQVW 


0.016 ± 0.004 


0.011 ± 0.003 


0.012 ± 0.007 


0.020 ± 0.016 


0.065 ± 0.035 


0.097 ± 0.064 


0.12 ± 0.22 


QV" 


0.005 ± 0.009 


0.018 ± 0.007 












VW" 


0.013 ± 0.011 


0.001 ± 0.008 












QVW" 


0.012 ± 0.007 


0.006 ± 0.005 












KaQVW" 


0.005 ± 0.005 


0.020 ± 0.004 













Note. — All entries have units of (fiK) 2 . The top set is for combinations of the pre-cleaned data. Sample variance is not included. The bottom 
set is for data cleaned with the KD3Pol model. Note that the cleaning has little affect on the 51 < £ < 150 bin other than to increase the 
uncertainty. 
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Table 9. Optical Depth vs. Data Selection 



Combination Exact EE Only Exact EE & TE Simple tau EE Simple tau, no I = 5, 7 



KaQV 


0.111 ±0.022 


0.111 ±0.022 






Q 


0.100 ±0.044 


0.080 ± 0.043 


0.08 ±0.03 


0.085 ±0.03 


QV 


0.100 ±0.029 


0.092 ±0.029 


0.110 ±0.027 


08S+ - 045 
u.uoo_ 015 


QV+VV 






0.145 ±0.03 


14+ 002 


V 


0.092 ±0.048 


0.095 ± 0.043 


09 +0 - 03 


1 0+ 003 
u - iu -0.07 


QVW 


0.109 ±0.022 


0.099 ±0.023 


0.090 ±0.012 


0.090 ±0.015 


KaQVW 


0.107 ±0.019 


0.105 ±0.019 


0.095 ±0.015 


0.095 ±0.015 



Note. - The values of simple tau are computed for 2 < £ < 11. The models are com- 
puted in steps of At = 0.005 and linearly interpolated. The last column is computed with 
the errors on I = 5, 7 multiplied by ten. The QV±VV is the QV combination without 
the QQ component. Since the exact likelihood is based on the Ka, Q, V, and W maps, 
there is no corresponding entry for QV±VV. Note that the maximum likelihood values are 
independent of frequency combination indicating that foreground emission is not biasing 
the determination of r. The calculations for the first two columns include the effects of 
marginalizatio n over synchrotron foreground emission and projecting out the small loss im- 
balance signal (I Jarosik et al.ll2006l ). 
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Fig. 1. — A model of the ionization history of the universe. The line marked "x" is the ionization 
fraction, x = n e /n where n e is the number of electrons and n = 11.2u;ft(l + z) 3 m -3 is the number 
of protons with w& the baryon density. Prom quasar absorption systems we know the universe has 
been fully ionized since at least z ~ 6. Between 6 < z < 30 the first g eneration of stars ionized the 
universe. We show a possible model inspired by [Holder et al.l (|2003l ). The history for this period 
is uncertain though the reionization produces a characteri stic signature i n the CMB polarization. 
For 30 < z < 2000, we show decoupling as described in iPeeblea ([1993 ). The line marked r is 
the net optical depth, t(z). The dashed curves are the integrands in the numerator (bottom) and 
denominator (top) of equation Q] (divided by 200) for the 100 < z < 2000 region. By eye, one can 
see that the ratio of the integrals at the maximum, and thus the fractional polarization, is ~ 5%. 
The vertical line marks the redshift of decoupling, z^ ec = 1088, at the maximum of the visibility 
function (not shown). 



-79- 



Tau A 



o 

CD 



CD 
CD 



CD 

Q 



1 2 5 10 20 5090 




187 



185 



183 



187 



185 



183 



187 



185 



183 



187 



Degrees (GC) 



Cen A 



24 



22 - 



O 
CD 



s 20 



CD 

Q 



18 



16 




185 



T(mK) 



183 



I ■ I ' I ' I ' I ■ 1 1 i I i I i I i I ■ I 



\ — / 1 / \ \ i 

- /♦■W-J*-- 

N /•/ **^-\ S 



^(J / 

. \ - / *r? \ \ 

. \ \ \ \ \ \ i. \ 
. \ \ \ \ \ \3to' 
' • / / / i fry' 

•*\\'\ U Ti~ 

,\\ \ ^ /,vy \* 
'^*< \-i IV 
.\»> J / /; 
I I 1 1 \ \ 1 I 
.Ml I \l^,\.f ^ 



; \ \ - 

v / «\ 1 
*^.\ > 

i < 

trrh, 
- I>*\^ 

\ > v~ 
— tj,] 

V V / - 
t\ .1 \ - 

1 \~ ' 

r-t'\ 1 
\ , 



312 310 308 312 310 308 312 310 

Degrees (GC) 



308 



312 



310 



308 



Fig. 1.— Top: Map of Tau A in Galactic coordinates at 41 GHz in Stokes /, Q, U, P, 
smoothed to 1°. Since Tau A is polarized parallel to the Galactic plane it is negative in Q 
and small in U. Bottom: Map of Centaurus A in Stokes /, Q, U, and P. For both sets of 
plots, Stokes I is scaled logarithmically and all the others are scaled linearly. The scaling in 
mK is indicated above the grayscale wedge for each panel. A map of the noise bias has been 
subtracted from the P images. 
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Q Band 



Fig. 3. — P and 7 maps for K, Ka, and Q bands in Galactic coordinates. See iBennett et al. 



(l2003bl . Figure 4) for features and coordinates. There is only one polarization map for K 
and Ka bands. For Q band, there are two maps which have been coadded. The maps are 
smoothed to 2°. The polarization vectors are plotted whenever a r4 HEALPix pixel (see 
§4.21 roughly 4degx4deg) and three of its neighbors has a signal to noise (P/N) greater 
than unity. The length of the arrow is logarithmically dependent on the magnitude of P. 
Note that P is positive. Maps of the noise bias have been subtracted in these images. 
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Fig. 4. — Similar to Figure [3] but for V and W bands. The two V-band maps have been 
coadded as have the four W-band maps. The relatively higher noise in the ecliptic plane is 
evident. Maps of the noise bias have been subtracted in these images. 
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Fig. 5. — A Lambert azimuthal equal area projection of the Galactic poles (left: north) 
showing the K-band polarization. The circumference of each map is at zero Galactic latitude. 
The convention in this plot is to use bars to indicate the polarization direction. It is clear that 
the polarization extends to high Galactic latitudes. A map of the noise bias is subtracted 
from this image. 
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Fig. 6. — Stokes Q and U maps in K and Ka bands. The Galactic plane is dominated 
by positive Stokes Q because the foreground polarization direction is perpendicular to the 
plane. As discussed in £JU this is expected because the Galactic magnetic field is predomi- 
nantly parallel to the plane. For comparison, the Stokes Q and U maps of a noiseless CMB 
simulation have peak-to-peak values of less than 6 fiK. These maps have been smoothed to 
1°. 
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o° 



180° 



Fig. 7. — Left: The angle of the magnetic field, 7m = 7PA + 90?, derived from the synchrotron 
radiation in the K-band map (smoothed with a 4° beam) shown in Figure 3. (We do not 
distinguish between ±180° in the field direction.) The predominant low Galactic latitude 
magnetic field direction is parallel to the Galactic plane (7m = 90°) and thus the synchrotron 
(and dust) polarization directions have 7 ~ 0°. In the North Polar Spur region, the magnetic 
field is perpendicular to the Galactic plane corresponding to 7m ~ 0° or 180°. Note the large 
scale coherency of the field. Right: The predicted magnetic field direction given by a simple 
model of the electron distribution and the logarithmic spiral arm model (Equation [9]) for the 
magnetic field. 



Reduction Factor for Synchrotron 



Reduction Factor for Dust 




Fig. 8. — The left panel shows the geometric suppression factor, g S ync{n), in the polarization 
due to the magnetic field geometry. The right panel shows a similar geometric suppression 
factor for polarized dust emission, gdustin), see §4.1 .31 
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Fig. 9. — Left: The observed K-band polarization, P. The color scale ranges from to 0.1 
mK. Right: The model prediction of the K-band polarization based on the Haslam intensity 
map. The model has one effective free parameter, the ratio of the homogeneous field strength 
to the total field strength as shown in Equation[13l This plot shows the results for f3 s = —2.7 
& q = 0.7. 
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Fig. 10. — Top: The H aslam 408 MHz map is shown with circles indicating loops from 



Berkhuijsen et al.l ( 1l97ll ). These ridges of enhanced Galactic radio emission are seen across 
the sky at low radio frequencies. The North Polar Spur ("Loop I") and the Cetus arc 
("Loop II") are examples of these features, which have been described as the remnants of 
individual supernovae, or of correlated supernovae outbursts that produce blowouts, or as 
helical patterns that follow the local magnetic fields projecting out of the plane. Four such 
loops can be seen in the Haslam 408 MHz radio map and the WMAP map. Note that 
the color stretch is logarithmic in temperature. Bottom: The WMAP K-band polarization 
map with the same loops superimposed. Note that the highly polarized southern feature 
is close to the North Polar Spur circle and may be related to the same physical structure. 
Note also that the polarization direction is perpendicular to the main ridge arc of the North 
Polar Spur, indicating a tangential magnetic field. This is also seen in the southern feature. 
Whether or not they are physically related remains unclear. 
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Fig. 11. — A map of the correlation, Z, between the polarization angle derived from the 
polarization of starlight, and the polarization angle in K-band. In the regions of high K- 
band polarization, the correlation is strong. The polarization directions are anti-correlated 
in the Orion-Eridanus region near I = —165°, suggesting spatially distinguished regions of 
dust and synchrotron emission. 
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Fig. 12. — The upper panels show the polarization signal at W band with the CMB and 
synchrotron signal removed (smoothed with a 10° Gaussian beam). The left and right panels 
show Stokes Q and U polarization components respectively. There is a clear preponderance 
of Stokes Q emission in the plane. The lower panels show the predicted dust polarization 
based on Equation [151 For \b\ < 10°, the stars do not sample the dust column well and the 
model is not accurate, especially for Stokes U. For |6| > 10°, there are regions where the 
data and model agree to the eye. However, a fit ( §4.31) is used to assess the level of polarized 
dust emission in the maps. 
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Fig. 13. — Temperature maps of the low polarization components for K, Q, and W bands. 
The maps are computed using equation (1161) . The color scale is in mK. Near the Galactic 
center, the low polarization component is approximately 6%, 3%, and 6% of the unpolarized 
emission in K, Q, and W bands respectively. 
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p02 p04 p06 p08 p10 Dust Sources 

Fig. 14. — The polarization masks, in Galactic coordinates, are shown for the P02, P04, 
P06, and P10 cut levels. The cross hatched region along the Galactic plane, common to all 
polarization masks, shows the dust intensity cut. The P06 cut is outlined by the black curve. 
The masked sources are in violet. The North Ecliptic Pole (NEP), and South Ecliptic Pole 
(SEP), and Galactic Center (GC) are indicated. 
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Fig. 15. — The Ka, Q, V, and W band Q Stokes Parameter maps before and after foreground 
subtraction using the method outlined in §4.31 There is a possible residual signal in W band 
though the noise is not yet sufficiently low to be certain. The U maps look similar. The 



clean ing for the cosmological analysis was done outside the processing cut (IJarosik et al. 



20061 ) and was based on the K-band maps and the starlight-based dust template. The over- 
subtracted dark regions on the galactic plane are inside the processing cut. 
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Fig. 16. — A comparison of the predicted Cg errors with (black) and without (red) assuming 
correlated noise in the polarization sky maps. On the y-axis is plotted the diagonal element 
of the inverse of the Fisher matrix for one year of data. The units are (/uK) 4 . Note that the 
y-axis scale for each plot is different. In each panel EE and BB are shown. The variations 
in the N _1 weighting are due to the scan pattern combined with the sky cut. There is 
less variation for B-modes than there is for E-modes. W4 has the largest 1/f noise of all 
radiometers. One can see that the combination of 1/f noise coupled with WMAP's scan 
strategy leads to a larger uncertainty than one would get from considering just the effects of 
1/f noise alone. 
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Fig. 17. — The absolute value of the EE (solid, violet through green) and BB (dashed, 
violet through green) polarization spectra for the region outside the P06 mask. The best fit 
ACDM model to TT, TE, and EE data with r = 0.09 and an additional tensor contribution 
with r = 0.3 is shown in black. The cross spectra have been combined into frequency bins 
according to Tableland into the following I bins: [2, 3, 4-5, 6-8, 9-15, 16-32, 33-101, 102-251, 
252-502]. In the presence of a dominant synchrotron spectrum, the averages over frequency 
are dominated by contributions from the lowest frequencies as can be seen by comparing the 
above at I = 2 to Figure [181 Diamonds (EE) and boxes (BB) denote the data points that 
are negative. The points are plotted at their absolute value to limit clutter. They should 
be interpreted as indicating the approximate noise level of the measurement. The la upper 
bounds and downward arrows mark points that are positive but consistent with zero. The 
general rise in the data for I > 100 is due to the large noise term. The red line corresponds 
to Equation [25] evaluated for v = 60 GHz for the BB foreground emission. 
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Fig. 18. — The frequency spectrum of the EE and BB power spectra for the region outside 
the P06 mask. To increase the signal to noise ratio, multiple values of I are averaged as 
indicated. Only statistical errors are shown. Negative values are not plotted. The frequency 
band combinations are given in Table The thin red line running close to the 1 = 2 EE 
spectrum is the model in Equation [261 The dot-dash red line corresponds to Equation [25] 
evaluated for BB at I = 2. 
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Fig. 19. — The frequency spectrum of the foreground cleaned EE and BB power spectra 
outside the P06 mask for £ = 2 — 9. Black shows EE, blue shows BB, and green, cyan, and 
orange show the EE yi — yj spectra (the BB ones are similar). For cosmological analysis, 
only the QQ, QV, and VV frequency channels are used (the "QV combination," indicated 
by red triangles on the bottom of each panel). The dotted black line shows the EE signal 
for r = 0.09. The dashed brown line shows the MEM dust temperature spectrum scaled 
by 0.0025 to indicate the level of 5% polarization (most evident near 90 GHz at I = 2). 
Averaged over the region outside the P06 mask, this is most likely an overestimate. The red 
curve shows the synchrotron spectrum scaled to 0.15 the pre-cleaned K-band temperature 
value. Based on the foreground model and discussion in text, it is unlikely that there 
is a significant residual foreground contamination in Q and V bands. Note that for all 
frequency combinations above 40 GHz (excluding KW), BB is clearly consistent with zero, 
also indicating the efficacy of the foreground cleaning. 
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Fig. 20. — Plots of all the noise for the expected null combinations of TB, EB, BB for 
the region outside the P06 mask. For T the foreground-cleaned V and W bands have been 
combined bands. For E and B, the foreground-cleaned Q and V bands have been combined. 
Cosmic variance is included for all plots. For each plot there are 999 I values that have 
been averaged into 33 bins for TB and 12 bins for EB and BB. For TB, x 2 / u = 41.6/33 
and x 2 / u = 931/999 with corresponding PTE = 0.15 and 0.94 for the two binnings. For 
EB, x 2 /" = 7.5/12 and x 2 / v = 956/999 with corresponding PTE = 0.82 and 0.84. For 
BB, x 2 /" = 6.2/12 and x l v = 1000/999 with corresponding PTE = 0.91 and 0.49. The 
polarization maps have been cleaned as described in Section §4.31 See also Table El 




Fig. 21. — Top: The EE and BB power spectra outside the P06 mask before and after apply- 
ing the KD3Pol foreground model. Different colors show different frequency combinations. 
Negative values are possible due to anticorrelations between foreground components, and to 
a lesser degree, from the coupling between different values of I. Only statistical uncertainties 
are shown. For EE, the smooth black lines are the best fit model to the TT, TE, and EE 
data. The cosmic variance uncertainty is indicated by the dashed lines. The EE values at 
t = 2 are 5.8, 4.5, & 5.5 jiK 2 for f/,40 (no KW), QVW, and QV combinations respectively. 
To clean these to a level of 0.1 /iK 2 requires cleaning the Stokes Q and U maps to one part 
in eight. The BB foreground emission is generally less than half the EE emission. Bottom: 
Expanded plots of the QV data for the P06 cut. The models are for r = 0.09 and r = 0.3. 
The absence of any signal in BB is another indication that foreground emission is not a 
significant contaminant. Note that the y-axis of the bottom plot has been divided by one 
power of I relative to the top plot. 
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Fig. 22. — The EE spectrum at I > 40 for all measurements of the CMB polarization. The 
curve is the best fit EE spectrum. Note that the y axis has only one power of I. The black 
boxes are the WMAP data; the triangles are the BOOMERanG data; the squares are the 
DASI data; the diamonds are the CBI data; and the asterisk is the CAPMAP data. The 
WMAP data are the QVW combination. For the first point, the cleaned value is used. For 
other values, the raw values are used. The data are given in Table [8] 
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Fig. 23.— Left: Lines of A% 2 for the fit given in Equation [28] vs. A EE . The different 
colors correspond to different frequency combinations. If the EE prediction from the TT 
measurements describes the EE measurements, then the minimum would be at A EE = 1. 
The line at A% 2 = 1 corresponds to the la error. One can see that the £ > 50 EE data are 
consistent with the model. Right: The amplitude and phase of the TE measurement with 
respect to the model predicted by the TT data. If the TE were completely predicted by the 
model based on TT, the contours would be consistent with A TE , A£ TE = (1,0). It is clear 
that the TT model describes the TE data as well. The reduced x 2 f° r the best fit model is 
0.67. To convert A£ to a phase angle in degrees, multiply by 1.18. 
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Fig. 24. — The TE power spectra for high and low I ranges for the region outside the P06 
mask in EE and TT. Left: At low I we use the QV combination for polarization with full 
N _1 weighting, and for temperature we use V band with uniform weighting. The black data 
points correspond to spectrum made with the KD3Pol cleaned polarization maps; the blue 
correspond to the same spectrum but witho ut cleaning (the I = 2 point is at 17.8 ± 3.4 /iK), 
and the brown are from lKogut et al.l (120031 ). The black dashed line is the best fit model to 
all the WMAP data. For the first-year data, x 2 = 35.3 for £ = 2 — 10 with a corresponding 
PTE « 0. For the three- year data, x 2 = 9-4 for £ = 2 — 10 evaluated relative to a null signal. 
The corresponding PTE is 0.4. When the three-year data are evaluated with respect to the 
best-fit model, x 2 = 5-4 with a corresponding PTE = 0.79. We find that the data sets are 
consistent with each other and that the three-year data prefer the r = 0.09 model over the 
null signal at the 2a (Ax 2 = 4) level. However, the three-year data are also consistent with a 
null signal. Right: The black data points show the three-year TE spectrum. This was made 
using V band for temper ature and the QV c ombination for the E-mode of polarization. The 
blue data points are from lKogut et al.l (120031 ). The smooth dashed curve is the best fit model 
to the WMAP data. An additional zero crossing near i = 400 is now present. 
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Fig. 25. — Plots of signal for TT (black), TE (red), EE (green) for the best fit model. The 
dashed line for TE indicates areas of anticorrelation. The cosmic variance is shown as a 
light swath around each model. It is binned in £ in the same way as the data. Thus, its 
variations reflect transitions between £ bin sizes. All error bars include the signal times noise 
term. The £ at which each point is plotted is found from the weighted mean of the data 
comprising the bin. This is most conspicuous for EE where the data are divided into bins 
of 2 < £ < 5, 6 < £ < 49, 50 < £ < 199, and 200 < £ < 799. The lowest £ point shows 
the cleaned QV data, the next shows the cleaned QVW data, and the last two show the 
pre-cleaned QVW data. There is possibly residual foreground contamination in the second 
point because our model is not so effective in this range as discussed in the text. The level of 
foreground contamination in rightmost two EE points could be roughly cr/2. For BB (blue 
dots), we show a model with r = 0.3. It is dotted to indicate that at this time WMAP only 
limits the signal. We show the la limit of 0.17 /iK for the weighted average of £ = 2 — 10. 
The BB lensing signal is shown as a blue dashed line. The foreground model (Equation 1251) 
for synchrotron plus dust emission is shown as straight dashed lines with green for EE and 
blue for BB. Both are evaluated at v = 65 GHz. Recall that this is an average level and 
does not emphasize the £s where the emission is low. 
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Fig. 26. — The relative likelihoods of r from the stand alone exact likelihood code and 
the first-year analysis. For the three-year results, all parameters except r and the scalar 
normalization, A, were held fixed as described in the text. The solid curves (labeled "WMAP 
3- years" ) show the exact likelihood for the Q V and KaQ VW combinations for the combined 
EE & TE data. If the K-band directions had been used for the dust polarization template 
(§4.3), leading to inferior cleaning, the likelihood curve would peak where the KaQVW 
does and have the width of the QV curve. The similarity indicates tha t any foreground 
contamination is small. The two broadest curves are from ISpergel et al.l (120031 ) and show 
the first-year likelihood for the WMAP data alone and for WMAP in combination with 
other data s e ts. T he dotted line is r likelihood for the first-year TE data as reported in 
Kogut et al.l (120031 ). The curve has a mean of r = 0.17 and width o = 0.04. 
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Fig. 27. — The two dimensional likelihood as a function of f and f for the BB spectrum. 
The contours indicate la and 2a. The n s parameter, which is degenerate with r and r, has 
been set at n s = 0.96. For the lightest contours, the tensor contribution to TT, TE, and EE 
is ignored. Because r is fully degenerate with r when the data are restricted to just BB, the 
limit is poor. The orange contours show the result when the TE and EE contributions are 
included, breaking the r — r degeneracy. The blui sh contours show th e result of including 
all data. The limit on r is more restrictive than in ISpergel et al.l (120061 ) because n s is fixed. 
When we marginalize over f , the 95% upper limits on f are 4.5, 2.2, and 0.27 for the three 
cases respectively. The plot shows that WMAP's ability to constrain r does not yet come 
from the BB data. The plot also shows that WMAP's ability to limit r depends critically 
on r. 
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Fig. 28. — The correlation coefficients of the Fisher matrices. The diamonds are derived 
from 100,000 Monte Carlo simulations, while the solid lines are the analytical formulae in 
the noise-dominated regime. In the simulations, the B-mode is noise only while the E-mode 
has some signals at low / from reionization. 



